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INTRODUCTION 


_ BE INVITED to deliver your Wright Brothers 
Lecture must at any time be reason for keen satis- 
faction, but in normal times would have set me the 
awkward question—am I competent to undertake the 
task? This year, when personal reasons tend to weight 
the scales, I confess that I have felt no hesitation—my 
desire to come was stronger than my doubts. You 
have done me great honor, and you have given me an 
opportunity, which every worker in aeronautics must 
covet, of paying my humble tribute to the authors of 
that great achievement now 38 years gone by. For 
both I am deeply grateful. 

It has become customary, in recent years, to preface 
aeronautical lectures such as this with lachrymose 
ponderings, brief or extended, on the question whether 
it was a good gift or a bad that Wilbur and Orville 
Wright gave to the world in 1903. I, who could never 
see human flight as a thing unconditionally good, do 
not propose to follow that tradition. The airplane 
was not the first gift of science which mankind has 
been free to use or to misuse, and it will not be the last. 
All power is double-edged, and always (unless the in- 
ventive faculty atrophies) we shall have to face re- 
sponsibility for the use we make of power. Let us 
ponder that responsibility, for it is both great and in- 
escapable; but not in scientific lectures. There let us 
rather give thanks, as men of science, for the fascinat- 
ing problems which aeronautics has brought us in its 
train. 

A New MATHEMATICS 


I was permitted to choose my title, and indeed, only 
on those terms could I have undertaken to lecture; 


* Presented before the Institute of the Aeronautical Sciences 
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for a man must talk of that with which his mind is: 
filled, and I confess that in the last six years my in- 
terests have been no longer general, my mind has 
been focused steadily and sharply on one single aim. 
That aim is the development of a new mathematics, a 
mathematics better suited than what is orthodox to 
the needs of engineering and aeronautical research. 
As a line of study it has proved most happily appro- 
priate to wartime conditions, with their pressure of new 
and urgent problems; and it still keeps our small 
“team”’ fully and contentedly occupied. Not all of 
our results may be revealed, some are still covered by 
the veil of wartime secrecy, but what can be released 
may perhaps serve to attract your interest and even to 
persuade some research worker to try our methods for 
himself. 

What then is my complaint against methods which 
are current? Chiefly this, that they take no cognizance 
of the unavoidable uncertainty of physical data. The 
physicists have their “uncertainty principle’: we too, I 
submit, should have our own. A simple example will 
serve to illustrate my point. (See Fig. 1.) 

Here is a “redundant framework,” redundant even 
when we neglect altogether the stiffness of its joints as 
distinct from the stiffness of its members. Similar ex- 
amples, but simpler, are expected by every candidate 
for examination in structural engineering, since the 
methods available for their solution are a commonplace 
of structural theory. What I want now to emphasize 
is a “bogus” feature of such problems as set in examina- 
tions, a feature which orthodox methods of solution do 
nothing to discourage. Only a very impractical ex- 
aminer would give marks for purporting accuracy (of 
calculated stresses or deflexions) closer (say) than 1 
per cent; the candidate, that is to say, may use a slide- 
rule. But the bogus feature is this, that the specified 
loading—5 tons here, 10 tons there—is specified as 














though the examiner knew it exactly—had it by supra- 
natural revelation. 

In fact, loadings are never known exactly. Usually 
aerodynamic or inertial, they are as uncertain as any- 
thing in the data of a problem, more uncertain (norm- 
ally) than the scantlings of the members or the elastic 
properties of their materials. To them, if to anything 
in the data, ‘‘tolerances’”’ should be attached: 5 + 0.1 
tons here, 10 + 0.3 tons there. But this is rarely done, 
even in cases we confront as problems of research; be- 
cause the doing of it would have no consequences, in ortho- 
dox methods of attack. 


THE MARGIN OF UNCERTAINTY 


Suppose that we do insert these tolerances—giving 
to every load the “margin of uncertainty” which it al- 
ways has in practice: then, clearly, any solution which 
accounts for the loads within their margins of uncer- 
tainty has the same claim as any other to be termed 
“correct.” ‘‘Exact’” solution becomes a word of no 
meaning in a theoretical treatment, just as “exact” 
measurement has no meaning in experiment. Will not 
theory gain reality from this parallelism with the ex- 
perimental outlook? 

For some 6 years I have sought, with a zealous band 
of assistants, to explore the implications of this new 
outlook. Breaking away from orthodox procedure, 
we have focused attention, always, not on the wanted 
quantities (stresses or displacements) but on the data of 
each problem, on the loads. We assume these to have 
always some specified ‘‘margin of uncertainty,” and we 
“whittle away” at them—accounting for more and 
more—until every ‘‘residual’”’ of a datum quantity has 
been brought within that margin. At this point we 
stop, asserting that to “whittle” further would be not 
only waste of time, but meaningless. If you cannot tell 
me the loads exactly (and in practice you never can), 
then you have no reason to expect unique solutions. 
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We need not apologize—we may even boast—that 
ours is “‘mathematics with a fringe.” 

Taking this standpoint—which I hold to be thus de- 
fensible philosophically—we have found our problems 
not made harder, but greatly simplified. This whit- 
tling process, if at times laborious, is easy: it could be 
taught to any intelligent bank-clerk. For it comes to 
this merely, that we give him a row of numbers (e.g., 
loads) which have to be accounted for, and a table of 
‘standard operations” (Table 1) whereby he can whittle 
them down to within a specified margin of uncertainty. 
Each operation, usually, is more powerful than any 
other to whittle some particular load: we indicate this 
in the ‘operations table,” and when “whittling’’— 
speaking generally—we go always for the biggest 
“residual” (the biggest load not yet accounted for). 
Errors need not be cumulative, and at the finish we 
have a complete picture of the wanted configuration— 
displacements, stresses, everything. Given time, we 
can work to any margin of uncertainty, however fine, 
so long as this is not zero. The wider the margin (and 
this «s left for the designer, not the computer, to assess), 
the quicker our route to a solution. Table 1 shows a 
typical ‘‘Operations Table” and the ‘‘Relaxation Table”’ 
in which this is utilized. 

Before the whittling process can begin, of course, we 
must construct the ‘operations table’: that is our 
first job—and surprisingly easy. For consider what 
is involved (again in relation to the typical framework 
of Fig. 1). Various operations can be imagined, but 
the simplest—yet quite general—is the imposition of a 
chosen displacement on a particular joint, all other 
joints being held fixed. Here is a problem special to 
each particular structure and to each particular joint; 
but it too can be left to our bank-clerk, if we can give 
him a general formula for the effect of any displace- 
ment imposed at one end of a single member when the 
other end is held fixed (Fig. 2). 


(a) (6) 
¢ 


THE UNIT PROBLEM 


So we come (for frameworks) to a “unit problem”’ 
soluble in general terms, and to formulas applicable 
without restriction. They are not greatly complicated 
when we take account of rigidity in the joints; and in 
their application redundancy—that bug-bear of con- 
ventional methods—no longer multiplies toil. 

Pursuing these ideas in relation to continuous girders, 
we found our treatment to be one, it essence, with the 
Moment-Distribution Method of Prof. Hardy Cross. 
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NEW PATHWAYS IN AERONAUTICAL THEORY 79 
TABLE 1 
(Units: 1 ton weight; 1 foot.) 
Operation | Nature of (1) | (2) (3) | (4) (5) 
no. operation X4 Xp vx Y, Yo X- Zz; Yp 
l(a) | u4=1 | —3,980 0 312 0 | 768 | 1,280 2,700 | —1,080 
1 () = 0.2513 | —1,000 0 785 | 0 | 193 322 678 —271s 
2(a2) | ug=1 i o | —4,100 | 0 (| 720 | 0 3,380 720 —720 
2 (b) = 0.2439 | 0 | —1,000 | 0 | 175 0 | 825 175 —175 
3 (a) 7% =1 | 312 | 0 | —5,393 4,500 461 768 — 1,080 432 
3 (b) = 0.1854 58 | 0 | 1,000 834 85s 142, | 200, 80s 
4 (a) vg = 1 | 0 | 720 | 4,500 —5,220 0 | 0 —720 720 
4 (b) = 0.1916 | 0 | 138 862 — 1,000 o | 0 —138 138 
5 (a) % = 1 768 | 0 | 460.8 | 0 | —2,480.8 | —768 0 2,020 
5 (b) = 0.4031 3095 | 0 186 | 0 | —1,000 —309; 0 814 
6 (a) %=%=1 | 312 | 720 — 893 —720 461 I 768 — 1,800 1,152 
6 (6) = 0.6200 | 193 | 446 —554 —446 286 476 ~—1,116 714 
7 (a) v4 = —0g = 1] 312 —720 — 9,893 9,720 461 768 — 360 — 288 
7(b) | = 0.0510 | 16 | —$7 —504 496 23 39 —18 —15 
Operation | Multiplier X 10° | X, | Xe | Yu | Yz Yo 
| (Initial forces) | | 15 | 7.07 | ae a 7. cs 
1 (b) 15 —15 | 0 1.18 0 2.89 
| ’ | Feel st te 2.89 
2 (b) 7 0 —-7.00 | 0O |} —1.22 0 
0 0.07 | 1.18 | 8.29 2.89 
6 (0) 10 1.93 4.46 | —5.54 | —4.46 2.86 
| 1.93 | 4.53 | -4.36 | 3.83 | 5.75 
7 (b) —8 | -0.13 | 0.30 | § 4.03 | -3.97 | —0.18 
1.80 | 4.83 —0.83 | -0.14 | 5.57 
5(b) | 7 | 2.18 | 0 |. 1.30 | Oo | —7.00 
3.98 | 4.83 | 0.97 | -0.14 | —1.43 
2(o) | 6 | oOo | -6.00 |} oO | 1.05 | 0 
|} 3.98 | -1.17 |. 0.97 | 40.91 | —1.43 
1 (6) | 5 | -5.00 | 0 0.39 | O | +0.97 
—1.02 | -—1.17 | 1.36 | 40.91 | —0.46 
6 (b) 2.5 0.48 | +1.12 | —1.39 | -—1.11 | 0.72 
| -0.54 | -—0.05 | -—0.03 | —0.20 0.26 
5 (6) 0.20 0.06 0 | 0.04 | O —0.20 
-0.48 | —0.05 | +0.01 | ~0.20 0.06 
1(o) | -0.5 0.50 o | -0.04 | 0 —0.10 
0.02 —0.05 | —0.03 —0.20 —0.04 
6 (6) —0.25 —0.05 | —0.11 | 40.14 | O11 —0.07 
/ -0.03 | -0.16 | +0.11 | —0.09 —0.11 
2 (b) | —0.2 | 0 | 020] 0 | -0.03, | 0 
| 0.03 0.04 | O.11 | 0.12, —0.11 
7 (0) 0.2 0 -~0.01 | -0.10 | 0.10 | 0 
| —0.03 0.03 | 0.01 | —0.02, | -0.11 
5 (6) | —0.12 ~0.04 o | -o02 | o 0.12 
| 0.07 0.03 | —0.01 | —0.02; | +0.01 
1() | —0.06 0.06 0 0 | oO —0.02 
| } 0.01 0.03 —0.01 | —0.02 —0.01 
6 (b) | —0.04 0 —0.02 0.02 | 0.020 | —0.01 
{| —0.01 0.01 0.01 | 0 | —0.02 





On simplifying assumptions which are customary, the 
same is true of rigidly-jointed structures. Here I 
would make full acknowledgment of Prof. Cross’ 
priority in the root idea of continued whittling or 
“liquidation”; and I also admit freely that my own 
notions (though they seemed to come out of a clear 
sky) may have had origin in his work, with which I 
had some acquaintance before I started, but no real 
familiarity. My notions came while I was writing, 
for a book, an account of methods for calculating 
stresses in braced frameworks: it was for that reason 


that I felt concern from the first to establish conver- 
gence, and it is to this circumstance that I attribute the 
ease of our path to harder problems, some of which I 
shall describe. There came to me suddenly, quite 
sharp and clear, a mental picture of a way in which real 
structures might be loaded, a way that could be followed 
in analysis. If I can convey that picture, I may leave 
it to you (I think) to convince yourselves that I have 
already given its mathematical equivalent. 

The ordinary ‘‘screw-jack”’ (e.g., for automobiles) is 
a means whereby a controlled displacement may be 
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TABLE 2 (N.B. All forces have been multiplied by 10°) (Continued on next page) 


AERONAUTICAL 





SCIENCES 


‘ail Line Imposed displacements d F, | F, F; F, | F; 
Stage 1 1 | Due to initial displacements (43) | 0.3 | —344,397 | 1,041,490 | —650,659 383,151 — 816,745 
2 | Due to displacement corrections (46) 376,829 | —501,980 672,426 | —350,339 549,192 


3 | Total forces | 32,432 | 539,510 21,767 32,812 — 267,553 
4 | a, (corrected) | 1.006674 | 1.596358 | 0.957343 | 0.449137 |—0.594696 
5 | Products (3) X (4) (Fy.a, and 7;.a,) 32,648 | 861,251 20,839 14,737 159,113 
6 | (F/T), = 0.08449; . > = 0.075887: 
F T | 
hence, according to (21), 9.21551 <A < 
0.224113 
Ad chosen for second stage = —0.077 0.223 
7 | Forces from line 3 (repeated) 32,432 539,510 | 21,767 32,812 — 267,553 
8 | Av.T, = —0.077T, —77,514 | —491,678 | —58,972 — 34,584 274,75 
Stage 2 9 | Corrected residuals (sums of lines 7 and 8) —45,082 | +47,832 — 37,205 —1,772 7,197 
10 | Displacement corrections (47) 37,595 | — 49,067 35,922 | — 279 — 2,928 
11 Total forces — 7,487 — 1,235 — 1,283 —2,051 4,269 
12 a, (corrected) 1.007154 | 1.627211 0.949575 | 0.440581 —0.602905 
13 | Products (11) XK (12) (F,.a, and T;,.a;,) —7,541 —2,010 | —1,218 —903 — 2,574 
14 | (F/T), = —0.000190; pe Pe —0,000971: 
F T 
hence, according to (21), 0.223190 <A< 
0.223971 
Ad chosen for third stage = 0.0008 0.2238 
15 | Forces from line 11 (repeated) —7,487 | —1,235 — 1,283 —2,051 4,269 
16 | AXN.T, = 0.0008T, 806 | 5,207 608 | 353 | —2,894 
Stage 3 17 | Corrected residuals (sums of lines 15 and 16) —6,681 3,972 —675 — 1,698 | 1,375 
18 | Displacement corrections (48) 6,495 —4,515 — 150 | 1,388 | — 641 
19 | Total forces | —186| —543} —825) ~—310 734 
20 | a, (corrected) | 1.004014 | 1.627319 | 0.949130 | 0.439763 |—0.602680 
22 | Products (19) X (20) (F,.a, and 7;,.a,) — 186 | — 884 | —783 — 136 | —442 


(F/T), = —0.0000835; PR Zz. =| 
ie 

—0.0001696; hence, according to (21), 

0.2238835 <» < 0.2239696 


imposed at some desired point: it is easy to imagine 
devices whereby this displacement can be recorded, 
together with the load sustained at any instant; also 
to visualize an arrangement in which every joint of a 
framework which would normally be free to move is 
provided with a jack of this kind to control its displace- 
ment. Suppose that initially the jacks fix the joints 
in positions such that the framework is not strained: 
then evidently, when the external loads are applied, 
they will be taken wholly by the jacks. Suppose that 
subsequently one jack is relaxed so that one joint is 
permitted to travel slowly through a specified distance: 
then load will be transferred from that jack to adjacent 
jacks and to the framework, and strain-energy will be 
stored in the latter. If the force on the jack which is 
relaxed had a component in the direction of the travel, 
that jack will be relieved of load, and strain-energy will 
be stored in the framework at the expense of the po- 
tential energy of the external forces. 

Here then is an imagined process whereby a frame- 
work could be brought from the unstrained to the fully 
strained configuration; and as I said, every step is easy 


to follow in calculation. It is a physical, indeed an 
engineering, picture of workmen instructed to “ease 
away” in turn, each man coming into the limelight as 
his recorded load becomes the biggest. 

All of this I have developed in a book published just 
a year ago, so I pass as quickly as possible to problems 
not yet written up (or at least, not in book form). 
But I would take leave for one exception, because with- 
out it I shall have difficulty in explaining more recent 
applications. That exception is vibrating systems 
fairly described, I think, as a problem of aeronautical 
interest. 


VIBRATING SYSTEMS 


I need not explain to engineers the significance in re- 
lation to crankshaft vibrations of problems such as 
Fig. 3 illustrates: we need at times to know the fre- 
quencies, and the associated ‘‘modes of vibration,” 
with which a given system can oscillate, left to itself. 
Nor shall I stop to elaborate ‘‘Rayleigh’s principle,” 
now widely known: it will merely give occasion for a 
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TABLE 2 (Continued) 
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F; ; Sums T; T2 T3 T, Ts T: Sums 
537,645 1,000,000 5,137,929 912,157 499,041 —2,517,769 —1,196,407 
— 576,392 6,674 1,247,503 | —146,283 —49,904 | —1,050,410 580,453 
— 38,747 1,006,674 6,385,432 765,874 449,137 — 3,568,178 — 615,954 
— 0.769943 1.006674 1.596358 | 0.957343 | 0.449137 —0.594696  —0.769943 
29,833 1,013,394 | 10,193,434 733,205 201,724 2,121,983 474,250 
(>> =1,118,422) (>> = 14,737,989) 
F 7 
—38,747 | 1,006,674 6,385,432 765,874 449,137 | —3,568,178 —615,954 
47,428 
8,681 | 
— 8,043 480 123,413 —6,214 — 8,556 — 49,251 1,470 
638 1,007,154 | 6,508,845 | 759,660 | 440,581 | —3,617,.429 —614,484 
—0.768106 1.007154 1.627211 | 0.949575 | 0.440581 | —0.602905 —0.768106 
—490 1,014,359 | 10,591,266 721,354 194,112 2,180,965 471,989 
~-* »~ 
oe = — 14,736 ) (> = 74,046 ) 
F T 
638 1,007,154 6,508,845 759,660 440,581 | —3,617,429 —614,484 
—492 
146 
37 — 3,140 429 — 356 —818 1,351 184 
183 1,004,014 6,509,274 759,304 439,763 | —3,616,078 — 614,300 
—0.767875 1.004014 1.627319 | 0.949130 | 0.439763 — 0.602680 —0_ 767875 
—140 1,008,045 | 10,592,663 720,678 193,391 2.179.336 71,705 


formula, the only formula that I shall need this even- 
ing. Here itis: 


V=p*T (p = 2xn) (1) 


We say that in “normal vibration”’ all parts of a sys- 
tem move in phase, so as to come to rest, and reach 
their maximum velocities, simultaneously: defines 
the frequency m, and then in Eq. (1) the left-hand side 
denotes the maximum potential energy (at the moment 
of rest), the right-hand side the maximum kinetic 
energy (at the moment of maximum velocity). Evi- 
dently these quantities are equal, since energy must be 
conserved: but now, suppose that we imagine con- 
straints to make the mode what we assert, suppose 
that we guess the displacements. Then Rayleigh’s 
principle, put briefly, tells us that we may be widely 
out in our guess about the mode, yet can use Eq. (1) to 
obtain a close estimate of the lowest natural frequency; 
and that our answer will never come out too low. Sup- 
pose then that by some other easy method we could 
arrive at a value known to be never too high: here 





would be knowledge almost if not quite as useful in 
practice as an “‘exact’’ result. 

Early last year we devised this second rule; and we 
saw, moreover, how to bring these vibration problems 
within range of Relaxation Methods, with their power 
of continuous adjustment. The result is a tabular 


computation (Table 2), 


involving as before operations 
performed according to the judgment of the computer, 
but in which, now, as the residual forces are brought 
within ever-narrowing margins of uncertainty, the 
upper and lower limits imposed on * continually ap- 
proach, like pincers holding the wanted value always 
between their jaws. The pincers never close, but we 
can make their gap insignificant, as the table shows. 
This is advance; but it would be advance of impor- 
tance only to the academic, if only the lowest frequency 
were required; for Rayleigh gives that, with all the 
accuracy that normally has value. The practical gain 
is this, that by our method we deiermine not only the 
lowest frequency but the corresponding mode; for thereby 
we are enabled to start afresh, this time to determine 


the second lowest frequency. In fact, for a system of 
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the kind that Fig. 3 exhibits, we can determine (and 
without excessive labor) al] of the wanted modes and 
frequencies. 

We have, moreover, found means of allowing for 
friction: we can, for example, deal with our former 
problem modified by the attachment of ‘vibration 
dampers” to two of the vibrating masses (Fig. 4). 
Table 3 gives (to those acquainted with these problems) 
the complete solution of this frictional problem. 


CONTINUOUS SYSTEMS 


Any method competent to deal with vibration prob- 
lems can also, of course, be applied to problems of 
elastic stability. We have found that ours are in 
fact thus applicable, but results may not be published 
now. I can say, however, that in all classes of problem 
we have found it possible to deal not only with systems 
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of finite freedom (such as I have described up to this 
point) but with continuous systems, such as are pre- 
sented in Elasticity and in Hydrodynamics. 

It is, in fact, in this field of two-dimensional problems 
that our work lies mainly now, and it is here that we 
are gaining our most encouraging results. We started 
in 1937, when we found it possible to solve (approxi- 
mately) not only plane-potential problems of an or- 
dinary kind (such as the torsion of non-circular bars— 
much studied of late by soap-film technique) but also 
problems which we had previously regarded as insol- 
uble—cases of plastic torsion, hollow shatts of non-cir- 
cular section, and the determination of magnetic lines 
in fields containing iron. Figs. 5-7 illustrate some 
published solutions of such problems. 

Fig. 5 shows our solution, by a synthesis of two particu- 
lar integrals, of the torsion problem for a pierced 
triangular shaft. The symmetry permits three dia- 
grams to be combined in one: the component solutions 
are exhibited at the bottom, and the final solution at 
the top of the figure. 

Fig. 6 exhibits the shear function for torsion of a 
triangular shaft as modified by the occurrence of 
“yield” at some limiting stress. Fig. 7 exhibits our 
calculation of magnetic lines in a field containing a 
triangular prism of iron, also (as a matter of historical 
interest) the solution obtained by Hele-Shaw and Hay 
(1900), with the aid of their very ingenious ‘‘stream- 
line-technique.” The correspondence will be seen to 
be very close. 

I should emphasize that in every case the triangular 
boundary was taken merely for illustration. What has 
been done in Figs. 5-7 could be done for any shape of 


boundary. 





TABLE 3 


Complete Solution for Small Damping of the Torsional Problem of Fig. 6 


y= (—0. 12583 B+i 0.16081 VA) \» = (—0.11126 B+i0.47325VA) | As = (—0.0090405 B+i 1.05471 VA) 
a; = 1+1 0.08692 B/ VA a, = 1.0040 =i 0.00476 B/ VA a4 = #£0.25832+7 0.07592 B/ VA 
a2 = 1.77927 +i 0.12229 B/ VA a, = 1.62732 i 0.07687 B/ VA a= £0.23515+7 0.06515 B/ VA 
a; = 2.56899+7 0.15196 B/ VA a; = 0.949137 0.02539 B/ VA a; = —0.83969+7 0.00809 B/ VA 
a, = 3.01088+7 0.11793 B/ VA ad; = 0.439767 0.03975 B/ VA a, = —1.036197 0.04004 B/ VA 
a; = 3.735387 0.08188 B/ VA a; = —0.60268+i 0.10722 B/ VA a; = —0.12846+7 0.02790 B/ VA 
ds = 3.83048+7 0.06864 B/ VA dg = —0.76787 +7 0.23549 B/ VA ad, = 1.890377 0.09874 B/ VA 

| 

\ = (—0.031401 B+ 1.14919 vy A) As = (—0.01239 B=: 1.56780 y A) | w= (—0.0017448 B +i 2.21649 vy A) 
a, = 0.26235+7 0.04002 B/ VA a4, = #1.34688*7 0.02737 B/ VA a4, = #£20.05797+71 0.00729 B/ VA 
a2 = £0.19510+7 0.01460 B/ VA a, = —0.22412+17 0.03730 B/ VA ad, = —0.123647 0.01589 B/ vA 
a3 = —0.91920 +7 0.07356 B/ VA a; = 0.01565+7 0.00050 B/ VA a; = 2.078967 0.00282 B/ VA 
a, = —1.00507 +7 0.09414 B/ VA 4,= 0.14098+7 0.02231 B/VA a, = —1.50198+7 0.00395 B/ VA 
a; = £0.29761+10.07491 B/ VA a; = —0.010347 0.00404 B/ VA ads = +0.04964+7 0.00778 B/ VA 
de = —1.11143+7 0.00782 B/ VA a¢ = 0.00761+:0.00276 B/ VA | de = —0.013367 0.00207 B/ VA 
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FINITE- DIFFERENCE APPROXIMATION 


These are all early attempts, and to explain them I 
must continue my account of methods. First, for the 
exact differential equation of our problem we sub- 
stitute an approximation in finite differences, thereby 
discarding the orthodox objective of determining a 
function everywhere, and setting ourselves instead to 
calculate its value at a large but finite number of chosen 
points, arranged symmetrically on a regular lattice or 
“net.” For our purpose each nodal point or ‘“‘knot” 
of the net must be surrounded by a fixed number of 
equal “strings” arranged symmetrically. There are 3 
possibilities, and only 3 (Fig. 8). 




















Fie. 8. 


Finite-difference approximations have been employed 
by various investigators—notably Sheppard and Rich- 
ardson in years preceding the first Great War. With- 
out exception, so far as I know, they have (in effect) 
entailed square nets (V = 4). We have found that 
triangular nets (V = 6) entail no additional difficulties 
and in general yield more accurate results. Hex- 
agonal nets (V = 3) do not appear to offer any advan- 
tage, but we have used the square and triangular indif- 
ferently, as each problem seemed to require. Some- 
times use of both is convenient, e.g., in our preparation 
(Fig. 9*) for the computation of flow of gas or steam 
through a convergent-divergent nozzle. Because we 
wanted to work on a square net, we were faced with 
the problem of conformally transforming the nozzle 
into a rectangle; but for that problem a triangular net 
was thought desirable on the score of accuracy. 

We have published (1940) a detailed account of our 
methods as applied to the problem of conformal trans- 
formation. There are four such transformations which 
give promise of being useful: a diagram will explain 
them better than many words. (See Fig. 10.) 


THE PLANE-POTENTIAL PROBLEM 


Though I have given examples from other fields, I 
need not labor, before this audience, the importance of 
the plane-potential problem to aeronautics. It is 
the first major problem in theoretical hydrodynamics, 
and must be solved if we would calculate stream-lines 
around an aerofoil or strut, must be solved as a prelimi- 
nary if we would calculate the effects of compressibility. 
Experimental devices, to obviate the need of calcula- 
tions, have been suggested frequently in recent years. 
It seems to have yielded unconditionally to computa- 
tion based on relaxation methods, and my young men 
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take it for granted, now, that any soJution they may 
need they can readily obtain for themselves, with all 
the accuracy that has value in physical applications. 
Lately we have been concerned to use our methods not 
as alternatives, but absolutely, for solving problems 
which by orthodox methods are intractable. I may 
not now give aeronautical examples, but perhaps it will 
suffice to give examples in the allied field of hydrody- 
namics. Here are three problems confronted in the 
theory of percolation—slow motion of fluid through 
the crevices of a porous material such as earth or sand. 


(Figs. 11-13.) 


THE ESSENCE OF THE RELAXATION METHOD 


But what, in these examples, is the part of relaxa- 
tion methods? I have told you only that we seek to 
compute the wanted function at points of a chosen 
“net”; I have not given details of the computations. 
Let us go back to first principles. 

Speaking quite generally, our problem is to deter- 
mine from its governing equations (in which I include 
boundary conditions) some wanted function z of two 
independent variables x and y. Here I would em- 
phasize that when I speak of a function I do not imply 
that its mathematical form is known: thus from my 
standpoint height of a country above sea-level is a 
function of position east-and-west (the x-direction) 
and north-and-south (the y-direction), but only a very 
courageous man would try to represent it in any detail 
by a mathematical expression. We exhibit its varia- 
tion in relief maps, or in contour maps which are a less 
elaborate alternative: similarly we could represent a 
mathematical expression, say z = xy, by a model in re- 
lief, but other means of representation exist. Here 
for example (Fig. 14) is a diagram from Jahnke and 
Emde which fulfils the same purpose exactly. 

Now, a new point. I have talked of models con- 
structed (or diagrams drawn) to represent known sur- 
faces or functional relations; but it is quite as natural, 
when the model is simple, to imagine it generated by 
And herein is seen the germ of 


some physical means. 
he first 


Prandtl’s well-known ‘‘membrane analogue’’: 


* War conditions made it impossible to send from England tlie 
original drawings from which the diagrams of this lecture have 
been abstracted, and the figures in consequence have had to be 
prepared from the author’s lantern slides. It is regretted that 
some have thereby lost definition of detail—notably Figs. 9-15. 
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pointed out that a membrane like a soap-bubble, un- 





der appropriate conditions in regard to pressure and the 








like, does in fact solve a plane-potential problem and 











et — will (as Taylor and Griffith were quick to demonstrate) 
Animale dremewsces se Pn te Ph 
r " solve plane-potential problems for us if we choose. 





I have said that in our computations we replace the 
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exact differential equation by a finite-difference ap- 
proximation, then seek values at the nodal points of a 
regular ‘‘net’’: now, we see physical reality in this 
hitherto mathematical device. The net of Fig. 8, pic- 
tured as strained by applied tensions, is our approxima- 
tion to Prandtl’s membrane; it satisfies the finite-dif- 
ference equation, and we can estimate the resulting 
error by iniuition. So in the end our way is clear: for 
we can treat the net just as at starting I treated the 
framework; we can (in imagination) furnish it like the 
framework with adjustable constraints; we can impose 
displacements on the nodal points, in any order we 
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Fic. 16. 


may choose; and throughout we may fix attention not on 
the wanted function (1.e., on displacements) but on loads 
not yet accounted for, seeking to bring these within 
chosen margins of error. That is the essence of the 
Relaxation Method. 

And in it lies, as I believe, the secret of our success: 
the use of relaxation patterns, typified in Fig. 15, to 
which are recorded no longer in 
Time alone 


dispose of ‘‘forces’ 
tables, but now on reticular diagrams. 
limits the completeness of their liquidation and the 


fineness of our net. 





Fic. 17. 
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What have we at the finish? No formal solution, 
merely a “grid” of computed values such as Figs. 
11-13 exemplify; coarse or fine as our effort has been 
short or sustained. This, let me emphasize, is what 
we must deduce from any formal solution before we 
can put it to practical use: the pure analyst, when 
having given us a formal solution he goes home to his 
tea, and from our standpoint, he has left his work half 
done! To any who object that it is not the whole pic- 
ture, I answer that it is what we accept without ques- 
ticn in a “half-tone’’ reproduction. Consider Fig. 16.* 

A. N. Black, a member of my staff at Oxford, has 
pursued this line of thought somewhat further, using a 
microscope to analyze part of a half-tone reproduction 
taken from a newspaper. Each elemental square was 
viewed in turn, and an estimate made and recorded of 
the fraction of its area which was “black”’ (7.e., inked): 
then, from the record of these estimates, a diagram 
was constructed on squared paper to have the same 
fractions inked of its corresponding squares. No 
particular convention was adopted—the inked frac- 
tion might have any shape or position in its square; 
but viewed from an appropriate distance the resulting 
diagram reproduces with surprising fidelity both the 
details and the general quality of the original photo- 
graph. Fig. 17 is, I think, self-explanatory. 


CONCLUSION 


I end, as I began, on a personal note. First, in 


* For the photographs here reproduced as Figs. 16 and 17, the 
author is indebted to Mr. V. Belfield. 


this lecture I have thought it invidious to name 
any of my research assistants, when I could not men- 
tion all: do not conclude that the hardest work has 
been done by me! I owe much to our Ministries of 
Aircraft Production, of Home Security, and of Supply, 
also to our Department of Scientific and Industrial 
Research, that they have given me the means to retain 
such zealous collaborators. 

Secondly, I could not even if I had wished have given 
this lecture a “war flavor,’’ for the very good reason 
that wartime applications may not be divulged; and 
it may well seem to you, in consequence, that through- 
out the last two years we have been working in aca- 
demic calm, seeking solutions for pure love of the quest, 
as is natural in times of peace. That is hardly a true 
picture of the facts, fortunate though we have been in 
the circumstances of our work at Oxford; for we have 
not let that work conflict with more immediate duties, 
and what I have given you tonight are but preliminaries 
and by-products of work which has throughout been 
focused on problems of the moment. What I have my- 
self contributed (projects for trial next day) has been 
written for the most part on my knee in railway 
trains—the only place, it would seem, where leisure can 
be found in England today. But I admit to having 
tried to act on the belief that isin me, that the course of 
new research cannot be controlled—the researcher must 
follow where his problems seem to lead. Let him try, 
with all the will that he can muster, to set himself the 
goal of usefulness to the national effort; but Jet him 
not (if he would obtain results) try too precisely to 
foresee the road. 


Discussion of the Lecture 


At the conclusion of the lecture, Dr. Hunsaker called 
upon the following for their comment and discussion: 


Dr. Hugh L. Dryden, National Bureau of Standards: 
“May I first offer my congratulations to the lecturer 
for the skill with which he has clothed the bare frame- 
work of a mathematical procedure of very general 
application, hidden its more forbidding aspects, and 
presented to us a physical structure with men hard at 
work with load-measuring devices and screw-jacks. 
The aeronautical world greatly appreciates Professor 
Southwell’s contribution in reducing the solution of 
many of its most difficult theoretical problems to 
addition, subtraction, multiplication, and division, 
processes with which most of us are familiar and which 
we can perform with or without the aid of computing 
machines. 

‘In this country we have yet to exploit, to the same 
extent as elsewhere, the fruitful union of mathematics 
and engineering in which each discipline adjusts itself 


to the limitations and needs of the other. The pure 
mathematician is inclined to regard the union as a kind 
of miscegenation, while the practical engineer’s relation- 
ship with mathematics is purely platonic. The greatest 
progress in aerodynamics, in structural theory, in 
meteorology, and in all other aspects of aeronautical 
science has come from mathematical theories guided by 
practical experimental results and from experiments 
interpreted in the light of mathematical theories. The 
empirical testing of a thousand specimens of an aircraft 
structure is worth less to the progress of aeronautics 
than one experiment on a single specimen carefully 
planned to determine which of the prevalent mathe- 
matical theories more nearly represents ‘the facts. 
For mathematics is able to lead us, as the lecturer has 
stated, along new pathways into the field of the un- 
known and untried—which is more easily and safely 
explored on paper than in the air. It is hoped that 
Professor Southwell will carry his gospel of mathe- 
matics with a fringe to the mathematicians as well as to 
the engineers. 
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“TI have no special qualification for discussing the 
more technical aspects which underlie the disarmingly 
readable account of the lecturer with its single mathe- 
matical formula in the text or to attempt to evaluate the 
relative convenience and merit of the relaxation method 
as compared with other methods of approximate nu- 
merical solutions of differential equations which are 
available. The lecturer has emphasized that it can be 
applied to problems for which the usual methods fail. 
Certainly one of the advantages of the relaxation 
method is that the bulk of the computation can be done 
by persons without special mathematical knowledge. 

‘Professor Southwell’s book gives a good picture of 
the services of the relaxation method in the theory of 
structures and in the theory of vibrating systems. 
Those of us with special interest in aerodynamics look 
forward to the time when aerodynamic applications, for 
example, to the flow in boundary layers and to the flow 
of compressible fluids, which cannot now be disclosed, 
are made available to us. 

“The members of the Institute are grateful to Pro- 
fessor Southwell for a most interesting lecture.”’ 


Dr. A. E. Lombard, Jr., Aircraft Branch, Office of 
Production Management: ‘‘Professor Southwell has 
given us a new approach to the general methods of 
successive approximation. Had it not been for the 
fact that he told us that this method was developed 
prior to the war, we would have suspected that he 
thought of it one night during a black-out in England. 
His method of relaxation sounds very much as though 
he had come out of a building on a dark night during a 
black-out and hit the first lamp post a rather violent 
bump whereupon he careened off, started down the 
street and hit the next lamp post a little less violently. 
Finally he has been able on this night to straighten out 
his course so that each successive lamp post down the 
street gets a lesser and lesser bump. 

“The British scientists have always been famous for 
their development of these methods of successive ap- 
proximation. Lord Rayleigh’s methods for vibration 
analysis and Hardy Cross’s methods for redundant 
structures are well known tous. Apparently, according 
to our lecturer tonight, the Southwell method is particu- 
larly suited in cases where there is an ample supply of 
bank clerks to carry out the work. The examples 
which he has shown us of the result of his method are 
stimulating and inspire us to investigate further the 
detailed workings of this particular method which shows 
definite suitability for the solving of certain problems. 

“Our lecturer also has emphasized a most important 
point to be realized by aircraft engineers and scientists, 
namely, that we are working with approximate data 
and our attempts at accuracy should be consistent with 
the quality of the data. There is an ever increasing 
attempt in the aircraft industry to carry out calcula- 
tions to a higher and higher degree of accuracy. This 
effort, when it is directed toward determining the 


fundamental physical constants with which we are 
working, is commendable, but when it is purely aca- 
demic juggling of arithmetic it is a deterrent to progress 
and may impede the design and production of new air- 


craft.” 


Dr. John E. Younger, University of Maryland: 
“T am truly astonished and amazed at three features of 
Professor Southwell’s presence among us: First, 
that he so nonchalantly makes the trip across the 
battle area of the Atlantic to honor us; second, that 
there is such close cooperation between his University 
and his government through him in this war-time 
emergency; and third, that he is succeeding in stream- 
lining structural analyses for war-time speed. The 
second and third of these features we will do well to use 
as our pattern. 

“Tf an engineer shows us how to build a new kind of 
airplane, a new kind of engine, or a new kind of wing 
flap, his contributions are always held with great ac- 
claim because everyone can see the material object and 
the benefit it will be to the science of airplane design. 
But when an eminent scientist presents us with the 
nucleus of a new mathematics much better suited to our 
aeronautical structural problems than our orthodox 
methods and the resulting development is obviously an 
enormous stride in the perfection of airplane design and 
structural analysis, unfortunately no such popular ac- 
claim may be expected. It is not spectacular, but it is 
the most needed improvement in the whole aircraft 
industry in this country. Structural design and struc- 
tural analyses are lagging far behind other phases of the 
industry. 

“‘There is too wide a gap between our mathematical 
elasticians and our structural engineers, so that the 
former, who are generally in charge of our research, are 
of little benefit to the latter, who are not concerned with 
ideal conditions and data but who must build cut of 
available stock, which, as Prof. Southwell has said, has 
an unavoidable uncertainty of physical data. 

“T have long felt that our methods of analysis of air- 
craft structures have been entirely inadequate for the 
purpose for which they were intended. They do not 
give the right answers. In general, the results of long 
mathematical analysis are no better for the practical 
structural analyses than an intelligent guess. We must 
bear in mind that even if an aeronautical structural 
engineer cannot solve the problem, he must neverthe- 
less find an answer which will be satisfactorily accurate. 
If we can obtain a satisfactory result for a structural 
design, taking all the factors into consideration, by a 
series of elementary approximations, it appears to me 
that we have made great progress in stream-lining our 
design stress analysis processes. 

“While we are on this subject, I wish to express the 
hope that Professor Southwell’s methods will enable us 
to satisfactorily explore the regions where Hooke’s law 


does not apply. Practically all our orthodox methods 


PRI ee 


ee ee 


Sages hegre? 


1 pate eer 





Re 


re at 


Aerie ab eS 





igo 


PEE SININIE EE OTRBR rT 


re 





‘ 
‘a 





PY 





NEW PATHWAYS IN 
are based upon Hooke’s law, but practically all our air- 
craft structural analyses are out of its field. The ortho- 
dox method of analyzing a structure within a propor- 
tional limit and assuming that the results apply above 
the proportional limit is always depressing because we 
know that our final result is nothing but a poor guess 
which may in general have been made without any 
previous analysis. 

‘‘We attempt too often to give our research reports a 
tone of scientific glamour by a long mathematical 
analysis which quite successfully conceals the value of 
the research from most engineers, but is not justified by 
the nature of the research or the results. It appears to 
me that Professor Southwell’s new pathways may 
quickly lead us to better coordination of our efforts.”’ 


Dr. George W. Lewis, N.A.C.A. (who was not pres- 
ent but who sent in the following comments): ‘‘The 
National Advisory Committee for Aeronautics ex- 
tends Professor Southwell a very cordial welcome to 
the United States, and appreciates that in these dis- 
turbed times it was possible for him to visit the United 
States and present the Wright Brothers Lecture. 
Members of the Committee’s staff are looking for- 
ward to the opportunity of having discussions with 
Professor Southwell, and it is hoped that his visit to 
the United States will be of sufficient length for him 
to have an opportunity to visit the Committee’s 
laboratories. ‘‘I regret that it was impossible for me to 
attend the Wright Brothers Lecture, but I did have the 
opportunity of reading the lecture. It has been said 
that the work of Hardy Cross on moment distribution 
is the greatest contribution to structural engineering in 
the present generation. The ‘new pathways’ of Pro- 
fessor Southwell hold the promise of being to engineers 
generally what moment distribution has been to the 
structural engineers. The important thing to me is 
that the methods developed by Professor Southwell 
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are not limited to aeronautics but can be applied in 
all branches of mathematical physics, of which engi- 
neering is one part. It is fortunate for us working in 
aeronautics to be among the first to have the benefits 
of this work. In the research of the National Ad- 
visory Committee for Aeronautics, relaxation methods 
have been used both in the solution of problems in the 
stability of structures and in the stress distribution 
of monocoque shells. The successful results obtained 
thus far indicate clearly that the use of this method 
will increase with time. It is regrettable that the 
present conditions do not permit Professor Southwell 
to give any details of his work on aeronautical subjects, 
which would, of course, be of particular interest to 
members of the Institute. Still, no matter how long 
we may have to wait for the release of this work, its 
intrinsic worth is such that the passage of time will not 
detract from its value.” 


Frank W. Caldwell, United Aircraft Corporation: 
“I take this opportunity to acknowledge the great 
debt which we owe to English science and particularly 
those branches of it that apply to aeronautics. Pro- 
fessor Southwell’s lecture brings to our attention the 
research done by himself and his colleagues which is 
of such fundamental importance to the progress of the 
aeronautical sciences. I want to thank Professor 
Southwell personally and as President of the Institute 
of the Aeronautical Sciences for the valuable contri- 
bution which he has made in the form of this lecture.” 

Dr. Hunsaker then presented to Professor Southwell 
the Wright Brothers Lecture honorarium of $500 which 
is provided by the Vernon Lynch Fund. He thanked 
Professor George B. Pegram of Columbia University for 
again making the facilities of the University available 
for the lecture. With an expression to Professor South- 
well of the appreciation of the Institute for his interest- 
ing and valuable lecture, the meeting was closed. 





On Some Vortex Theorems of 
Hydrodynamics 
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DYNAMICS OF INCOMPRESSIBLE FLUIDS 


HE MOTION OF each element of an incompressible 

fluid takes place in coordination with that of all 
others in such a manner that no element of space be- 
comes incompletely filled or overcrowded. The mathe- 
matical expression for this kinematic condition is the 
equation of continuity. 


V-v=0 (1) 


where v denotes the velocity of the fluid. This equation 
expresses the fact that the net outflow of fluid over any 
infinitely small closed surface, divided by the en- 
closed volume, is zero. This equation applies to every 
incompressible fluid either with or without viscosity, 
for steady and unsteady motions alike. 

The mechanism enforcing this continuity of flow is 
a system of internal stresses between adjacent fluid 
elements. These stresses are called into being as soon 
as any influence comes into play tending to compress 
or rarefy the fluid. For non-viscous or ideal fluids this 
system of stresses consists simply of a distribution of 
normal stresses or pressures. 

The basic problem of hydrodynamics is the deter- 
mination of the fluid motion when the boundaries of 
the fluid and the external forces acting upon the fluid 
are given. In many applications the converse problem 
is encountered, namely the determination of a system 
of external forces which is necessary to bring about a 
given motion of the fluid. It is important to realize 
that the latter problem does not have a unique solu- 
tion. There are an infinite number of possible force 
systems which result in the same fluid motion. In 
order to further limit the force systems it is necessary 
to impose additional restrictions regarding the nature 
of the force systems, such as, for instance, that the ex- 
ternal forces act only through specified regions. The 
reason for this is that certain force systems merely tend 
to compress or to rarefy the fluid without tending to 
The action of gravity on a glass 
Such ex- 


influence its motion. 
of water standing immovable is of that kind. 
ternal forces are counteracted by a modification of the 
pressure distribution; they are purely static and do 
not influence the fluid motion. 

As a first step in an analysis of fluid motions such 
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purely static force systems and their resultant static 
pressure distributions are dismissed from the investi- 
gation. Two force systems will be considered as equiva- 
lent if they have the same dynamic effect. Two equiva- 
lent force systems f, and f, will therefore differ by a 
purely static force system which can be balanced by a 
distribution of pressures P. This fact can be mathe- 
matically expressed as follows: 


f. —f,; = VP 


The negative gradient, called the anti-gradient of the 
pressure distribution P, represents the buoyancy force 
per unit volume incidental to the pressure distribution 
P and in this case balancing the static force fp — f). 

The rotation of a gradient is always equal to zero, and 
thus 


VxvP=0 
Therefore 
Vv X (& — f;) = 0 
or 


VXfh=V Xf, (2) 


The condition for the equivalence of two force systems 
is therefore that they have equal rotation (curl) at all 
points. This statement holds also for viscous fluids 
because the static stress system of an immovable vis- 
cous fluid is likewise a pure pressure distribution. 

This introduction is concluded by remarks on the 
external forces encountered in hydrodynamic prob- 
lems. The gravity forces acting on a fluid of constant 
density form a static system and can of course be taken 
care of by a suitable pressure distribution in conjunc- 
tion with suitable boundary conditions. With ordinary 
reference systems, there remain then as dynamic ex- 
ternal forces acting on the fluid, those transmitted by 
foreign objects immersed within the fluid or forming 
the fluid container. If the objects are infinitely thin 
sheets or lamina, their effect can be entirely described 
by disregarding the lamina, considering the fluid as 
continuous, and considering only the forces transmitted 
by the lamina to the fluid. If the immersed objects 
are voluminous, or constitute the external fluid bound- 
aries, they may likewise be disregarded, and the space 
occupied by them may be assumed to be filled with 
fluid at rest and of the same constant density as the 
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fluid investigated. The reaction between the solids 
and the fluid is again described by considering the 
respective forces on the fluid only. In that way, every 
flow problem can always be reduced to an equivalent 
problem referring to the entire space being filled with 
the fluid in question and acted on by a system of exter- 
nal forces. 


THE ACCELERATION 


Consider the motion of a perfect fluid (no viscosity). 
The motion may be steady or unsteady. Then the 
problem is to determine the force system, or more ex- 
actly, one of the many equivalent force systems which 
is required for the given fluid motion. According to 
Newton’s Second Law the force f acting upon a small 
element of the fluid having unit volume is equal to the 
product of the mass p and the acceleration a of the ele- 
ment, that isf = pa. The mass p of the unit volume is 
called the mass density. The acceleration a may also 
be expressed as 


a = Dv/DT = (dv/OT) + v-V; v (3) 
(See reference 2, page 4.) 


The second term in this expression may be transformed. 
From 


vX(V Xv) =v-V; v—vv; V 


=v-V¥V; v— '/Vv-v 
it follows 
a = Ov/OT + '/.Vv-v + Vv X (V X v) (4) 


From what has previously been discussed, the chief 
interest is in the rotation of the external forces, and 
hence in the rotation of the acceleration. This is equal 
to 


V Xa=0u0/0T+V7 X (v Xu) 


where u = V Xv. The term V X V (v-v) vanishes 
since the rotation of a gradient is zero. 

For steady flows the term 0u/O7 vanishes by defini- 
tion. The rotation of the acceleration then reduces to 


Vx<Xe=V XC Xp» 
and thus 
VxXf=pV X(vXv) (5) 


This leads at once to the basic theorem about the force 
distribution required to keep the motion of a perfect 
fluid steady. 





Basic Vortex Theorem 


A force distribution equal to the cross product of the 
velocity and of the vorticity, multiplied by the fluid density, 
maintains any motion of a perfect fluid as a steady motion. 

This is not Kutta’s theorem. The forces of airfoil 
theory are directly . presented by the basic vortex 
theorem. The forces transmitted by an ideal propeller 
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disc to the ideal fluid are not directly specified by the 
above theorem, but are merely equivalent to the 
specified forces. 

In the absence of external forces, the force system of 
the Basic Vortex Theorem must vanish. 


Steady Flow Criterion 


The condition that the initial motion of a perfect fluid 
be steady without the action of external forces is that the 
cross product of the vorticity and of the velocity be irrota- 
tional. 

This statement should be compared with the dif- 
ferent but not inconsistent statement in reference 2, 
page 226. 


HELMHOLTZ VORTEX THEOREM 


An important progress of the theory was made by 
considering the rate of change of the circulation of the 
velocity along a closed path moving with the fluid, 
such that the path remains occupied by the same fluid 
particles. This integral is a Lagrangian quantity and 
has to be expressed in Eulerian terms. The integral 
considered is 

D/DT f dr-v 
This is equal to 
D/DT f dr-v = { (Ddr/dT)-v + fdr-(Dv/DT) (6) 


the first right-hand term, being {“dv-v = '/2v-v/, is 
clearly zero, the integration being performed along a 
closed path and returning to the point from which it 
started, at which v-v must have the initial value. There 
remains therefore only the second right-hand term. It 
follows that the rate of change of the velocity circula- 
tion is equal to the circulation of the acceleration. 
That, multiplied by the density, gives the circulation 
of the forces. If these forces have no circulation, that 
is, if they are irrotational, the velocity circulation along 
a path moving with the fluid will not change. 


Helmholtz’ Theorem 


In a perfect fluid, not under the action of external forces 
having rotation, the vortices move with the fluid. 

Conversely, it requires such external forces to bring 
about fluid motions in which there is a relative velocity 
between the vortices and the fluid elements; and the 
larger this relative velocity is, the larger will be the 
forces required for equal vortex strength and density. 


GENERALIZATION OF THE BASIC VORTEX THEOREM 





The Basic Vortex Theorem refers to steady motions 
only. Unsteady motions may now be considered, in- 
quiring about a similar vortex theorem that may ease 
the burden of the mathematical investigator. 

The most general unsteady flow cannot be brought 
about by external forces acting at regions of vorticity 
only, and thus being simply associated with the vor- 
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tices. For the terms Ov/OT do not stand in any neces- 
sary or direct relation to the vortex system u. 

A direct relation between the forces and the vortices 
can, however, be pointed out for the special case that 
all vortices can be considered to move merely, with a 
velocity w, without the creation or annihilation of any 
vortices. Each point of each vortex line is transferred 
through the distance w, the new points are connected 
with each other to form the new vortex line, and the 
circulation or strength initially ascribed to the old 
vortex line is again ascribed to the new one. Let there 
be only one such isolated vortex line of concentrated 
vorticity. The forces can then easily be determined 
from the theory of potential flows. The strip composed 
of surface elements of width w and parallel to u X w 
extending along the vortex is then the only region where 
the potential of the difference of the initial and of the 
modified velocity distribution abruptly changes its 
value. This change is equal to the circulation of the 
vortex, u. Hence a force p(u X w) is necessary to bring 
about the change. The force for maintaining the steady 
flow was seen to be p(u X v). Hence the total force 
pu X (v — w) results. 


Generalized Vortex Theorem 


The unsteady motions of perfect fluids require external 
forces equal to the cross product of the vorticity and its 
velocity relative to the fluid, multiplied by the density. 

This very special demonstration is nevertheless 
conclusive of the general validity of the theorem. 

Proceeding to a general proof of the generalized vor- 
tex theorem, it is known that the rotation of the force 
responsible for the unsteadiness has to be the space 
integral of pOu/O7. If any vector is free of divergence, 
its space integral is equal to the surface integral of its 
flux multiplied by the radius vector, that is to say it is 
equal to the static moment of that flux. This follows 
from 

Vu; r—V-ur=uV;r=u (7) 


Sdo-ur= fdSu 


(See also Eq. (28) of reference 3.) The space integral 
of a system of vectors free of divergence can therefore 
be expressed from the surface conditions of the region 
considered, or, otherwise expressed, vector modifica- 
tions leaving the surface flux intact leave also the 
space integral intact. 

Consider a region of the space and substitute for the 
vortex configuration prevailing after the motion w 
of the vortices another, modified one, namely the initial 
vortex configuration present before the motion w took 
place. To this system are superposed near the surface 
two incomplete vortex systems: (a) the vortex ele- 
ments carried through the motion w from the outside 
into the inside of the region, and (b) vortex portions 
of original strength and direction, connecting the interior 
ends of the moved vortex portions with their original 
point of intersection with the surface. 
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This modified vortex system, free of divergence, 
gives the same integral as the vortex system actually 
prevailing after the motion w. Compute the space 
integrals of the two incomplete vortex systems (a) and 
(b). As to (a), the motion w carries into the space 
element considered, through its surface element do, the 
vorticity — (do-w)u, the expression within the bracket 
being the flux of the velocity w. The system (b), vor- 
tices of extension w and strength do-u, give elements 
(do-u) w. Hence, the total change of the vorticity 
will be the surface integral of 


do-u w — do-wu = do X (u X w) (8) 


Reducing the space element considered to infinitely 
small size, then carrying out the surface integration (8) 
and dividing the outcome by the volume of the space 
element, leads to V X (u X w) by the very definition 
of the rotation. This corresponds to a rotation of the 
unsteadiness force pV X (u X w). 

The rotation of the steadiness force was pV X 
(v X u) according to the basic vortex theorem. Hence 
the rotation of the total force is seen to be equivalent to 
the rotation of the product of density times the cross 
product of vorticity and the relative vortex velocity. 


APPLICATION OF THE GENERALIZED VORTEX THEOREM 


It has already been pointed out that the generalized 
vortex theorem relates to special flows only. It cannot 
be generally applied to each and every unsteady flow of 
a perfect fluid. The application is particularly called 
for if all points of each vortex line move with an equal 
relative velocity. This is always the case with two- 
dimensional unsteady flows. The theorem is then ap- 
plicable. It is, for instance, applicable if all such vor- 
tices occupy a straight line and continue to occupy the 
same, the velocities at the points of the straight line 
being parallel thereto, and hence the forces being at 
right angles thereto. These conditions prevail in the 
investigations reported.‘ 

In that case, the relative velocity of the vortices is as 
follows. Let x denote the Cartesian coordinate parallel 
to the line: 

ou/OT = w(du/dx) (9) 


and hence 
w = (0u/07)/(0u/dx) (10) 


Hence the force per unit length of the line, necessary to 

bring about the unsteady flow considered is equal to 
piv — [(0u/d7)/(du/dx)] XK u} (11) 

By the integration of this force and of its moment, the 


total force and the total moment are directly obtained. 


FURTHER GENERALIZATION OF THE BASIC VORTEX 


THEOREM 


The limitation of the generalized vortex theorem to 
certain modes of fluid motion can be removed and, 
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formally at least, the theorem can be applied to the 
most general case. So far, it is necessary that there are 
vortices present in the first place, which by their rela- 
tive motion may provide the force reaction. If vortices 
are absent, they may however be provided by the in- 
troduction of fictitious vorticity. This is a pair of sys- 
tems of equally strong but opposite vorticity, super- 
posed to each other and thereby canceling each other. 
The proposition is indeed academical, and probably 
devoid of technical usefulness. But it is highly instruc- 
tive and of heuristic value. A larger completeness and 
coherence of the science and its geometric pictures is 
brought about. It can thus be demonstrated that the 
generalized theorem leads back to the basic Eq. (5). 

Consider a space element dS at rest, bounded again 
by surface elements do. Let there act a force system f 
on the fluid occupying and surrounding the space ele- 
ment, said fluid moving there relative to the space 
element with the velocity v. Introduce now such fic- 
titious self-canceling double vorticity parallel to the 
surface elements and at right angles to the force. In 
accordance with the generalized vortex theorem, the 
two vortex systems are then moved across the surface 
elements, in opposite directions indeed, but also with 
opposite intensities, and therefore equally contributing 
to a transport of vorticity into the space element. 

Before computing this, disregard those components 
of the forces already required for the maintaining 
of the steady motion, namely p(v X u) according to the 
basic theorem. The force to be accounted for by the 
motion of the fictitious vortices is accordingly f — 
o(v Xu). This force gives rise to a transport of vorticity 
through the surface elements 


du/dT = fdo X (f/p — v X u) 


Reducing the space element to small dimensions and 
unit volume and carrying out the integration gives 
finally 


du/dT =V Xf/p -V X (VX 4) 


which corresponds exactly to Eq. (5), by substituting 
f/p for the acceleration. 


INTEGRAL THEOREMS. THE MOMENTUM 


The preceding theorems deal with the relation be- 
tween the forces at each individual point and the vor- 
ticity at the same point. The following theorems are 
of an altogether different class, dealing only with the 
total or integrated value, one single numerical value or 
set of values, either of the momentum, or of the moment 
of momentum, or of the kinetic energy of a body of a 
perfect fluid. The substitution of the vorticity for the 
velocity itself reduces the computation work in all 
cases where the vorticity is restricted to surfaces or 
lines. Thus a space integration may be reduced to a 
surface or line integration. 


First the momentum. The impulse applied to the 


fluid filling infinite space, on the occasion of bringing 
about a certain flow configuration, is not recognizable 
subsequently in the total momentum of the infinite 
body of fluid. During the creation of the motion, there 
occur in general impulse exchanges between portions of 
the fluid, disproportionally surpassing in magnitude the 
impulse applied externally, and rendering it unrecog- 
nizable. Portions of the fluid have acquired an infinite 
momentum in one direction, other portions an almost 
equal momentum in the opposite direction. The sum- 
ming up process leads to an ambiguous result, because 
this result depends on the order in which the contribu- 
tions of the various portions of the fluid are successively 
taken into account. 

While thus a momentum of the infinite body of fluid 
is not uniquely ascertainable and insofar does not 
exist, it is nevertheless of practical value to have avail- 
able one unique preferred value of the total momentum, 
for the computation of the impulse required for the 
creation of the motion, and for some computations re- 
lated thereto. A method how to define and obtain 
such preferred value is proposed and explained in refer- 
ence 3. The momentum is computed from the fictitious 
irrotational motion of similar fluid imagined to be in- 
closed within the surface of the immersed solid, which 
motion agrees with the actual outside irrotational 
motion at all points of the surface. Of course, for a 
finite body of fluid there is no inherent difficulty in de- 
termining the momentum. There is, however, a particu- 
larity of the fictitious flow not found in physically real 
flows. In spite of the assumed constant density, and 
not inconsistent therewith, the flow may possess di- 
vergence; the fluid must be supposed to appear or dis- 
appear at certain regions. 

The existing theorems: serving for the computation 
of the momentum and the related quantities in case of a 
finite fluid body are therefore in need of a revision and 
particularly exact statement, to establish whether they 
also hold for the case of divergence. The variety of 
theorems will at the same time be increased by a few 
similar theorems which are possibly new and relate to 
the moment of momentum. At present, however, it 
will be necessary to dispose of the momentum itself. 

The momentum of the fluid contained within a cer- 
tain region of space has. been expressed in reference 3, 
Eq. (28) by means of its divergence or source and sink 
system including its boundary flux system. It is there 
shown that the total momentum is equal to the static 
moment of the divergence and of the boundary flux 
multiplied by the density. 

This integral gives the total momentum regardless 
of the presence of or coexistence of any vortex system, 
including the ‘‘roll,’’ the surface integral of the cross 
product of the surface elements and the velocity. That 
is to say, the momentum is expressed in terms of the 
divergence only, for the boundary flux is equivalent to 
and complimentary to the divergence, becoming genu- 
ine divergence as soon as the space beyond the boundary 
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is supposed to be filled with resting fluid and included 
in the integration. 

The static moment of the divergence is the prototype 
for other expressions containing the rotation and valid 
nevertheless for any coexisting divergence distribu- 
tion. The disregard of either the rotation or the diver- 
gence in either case is only formal, and the disregarded 
quantity not explicitly contained in the integral ex- 
pression still influences the final result. For it is still 
associated with the boundary conditions, and so makes 
itself felt in the boundary term if the truth requires 
that. If the fluid fills the entire space, the external 
regions are supposed to be occupied by resting fluid, 
and this presupposes and requires certain mutual rela- 
tions between the divergence and the rotation or its 
absence by virtue of which certain integrals of one are 
also determined by the other. 

The entire class of integral theorems now under dis- 
cussion have this in common: that their surface 
terms stand in a unique relation to their space terms, 
in that they are fundamentally the same specially trans- 
formed for surface effects. The divergence is replaced 
by the flux, and the rotation is replaced by the roll. 
The surface terms appear from the space terms by par- 
tial integration or application of Gauss’ theorem, 
which constitutes a generalization of Green’s theorem. 

The possibility of these theorems is intimately con- 
nected with the fact that V designates a differentiation 
dimensionally equivalent to a division by a length. 
This division effect can be counterbalanced by a suit- 
able multiplication by r, the radius vector. 

All this is most clearly exhibited in the divergence 
integral for the momentum. The space term (V -v)r, 
is dimensionally correct, being a vector; and one fac- 
tor r counterbalancing the differentiation V. Only 
v is differentiated, r is not. The surface term do-v r 
results from the space term by replacing V by do, or 
the divergence by the flux. The significant part of the 
moment, the space integral over v, results from the 
above space term by dislodging V from v and connect- 
ing it in direct association with r, preserving otherwise 
the vector-analytical equivalence. 

The same sequence of steps may now be applied to 
the expression 


(Vv Xv) Xr 
which is a vector of equal geometric dimension as 
(V-v)r 


The complemental surface term is (do X v) X rf, 
which can be written down without any analytical 
effort. The turning over of the differentiation V from 
v to r is a little more complicated than before, as long 
as one is restricted to ordinary dyadics. It is then 
necessary to dissolve the expression first to make all 
cross products disappear. 


TM ¥) Xe 


V;vr-—-v; V-r 
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Its complement: 
V;rv—vV-r= —2v 
Hence 
SdS(V Xv) Xr—2fdSv= f{ (do Xv) Xr (12) 
SdSv = '/2fdS(V Xv) Xr —/2f (do X v) X r (13) 


Momentum Theorem 


The momentum of the fiuid occupying a region of space 
zs equal to the sum of half the vector moments of its vorticity 
and of its surface roll, multiplied by the density. 

The factor '/2 can easily be verified for the case of a 
very small cubical space element being occupied by fluid 
flowing parallel to one edge. The space integral be- 
comes then zero. The surface integral, without the 
factor '/2, gives then twice the volume multiplied by 
the velocity. 

The theorem is also set forth in reference 2, Section 
152, in a different form, not inconsistent with the one 
just given. The factor '/2 is there also correctly stated. 


THE MOMENT OF MOMENTUM 


No suitable expression containing the divergence 
only can be found for the expression of the moment of 
momentum. It is indeed apparent from elementary 
considerations that the same divergence may be as- 
sociated with fluid motions having different moments of 
momentum. . 

Three expressions present themselves as a first choice 
for obtaining integral expressions for the moment of 
momentum containing the rotation only. These are 


(a) V X vr-r, the moment of inertia of the rotation 
(b) ((V Xv) Xr) Xr the vectorial moment of in- 
ertia of the rotation, and 
r-V Xv)? 
(c) (WT XWXDXr-(v XW 


VV xXvVV¥ xX¥ 


an apparently more complicate expression, but adapt- 
able to easy expression in Cartesian coordinates. 


Case (a) 
SdSV Xvrr= fdoXvr-r+ fdSv X (Vr-r) 


= fdo Xvrr+2fdSv X r(14) 


Hence 


SdSvX r= V2 fdS(V Xv) r-r— 1/2 f'do X vr-r (15) 





First Moment of Momentum Theorem 


The moment of momentum of a body of fluid is equal to 
one-half the sum of the common moment of inertia of its 
vorticity and of the same moment of momentum of its sur- 
face roll, multiplied by the constant density. 


Case (b) 


(V Xv) Xr)Xr=(V Xvrr-—VXvrr (16) 
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JSdo Xvrr+2fdSv Xr —- 
SdS(V Xv)-rr (17) 
=—fdSrXv-V;r— fdoXvrr— fdSv-V Xrr 
= f[dSvxXr- fdoXvrr—o (18) 


SdSV Xvyrr= 


whence 
SdSv Xr = —fdSV Xv)-r; r+ fdoXvrjr 
(19) 

which expresses the moment of momentum in terms of 
the dyadic moment of inertia of rotation and roll. 
Combining (15), (17) and (18) 
SdSvXr='%3fdS((V Xv) Xr) xX rt+ 

/3f ((do Xv) Xx) Xr (20) 
Second Moment of Momentum Theorem 


The moment of momentum of a body of fluid is equal to 
one-third of the sum of the vectorial moment of inertia of 
its vorticity, and of the same moment of inertia of its sur- 
face roll, multiplied by the constant density. 
Case (c) 

The first term is the same as under (a). The second 
term, 


(r-V Xv)? 
Vv ————— 
V Xv: 


VX 


leads to a complement containing as a factor (r X V -v) 
and is accordingly zero. Hence the factor '/2 appears 
again. 

This third theorem is stated in a different form, but 
consistent with the form here given, in reference 2, 
Section 152. A numerical example for its application 
is there also given. No references are known to the 
author for the first two or three moment of momentum 
theorems. 


THE KINETIC ENERGY 


Flows free of divergence and flux, as well as flows free 
of rotation and roll, possess each a definite kinetic 
energy. It becomes therefore apparent that the ex- 
pression for the integral must contain both divergence 
and rotation. It must furthermore be of second order 
in either derivative and the velocity. 

It is sufficient for the application that each term con- 
tain the rotation or divergence one time as factor only, 
because this is enough to restrict the integration to 
limited regions if that can be done at all. 

An expression to start from, which leads to the aim is 
as follows: 

VXverXv4t+V-vr-v= (21) 

(Vi v-r)-v—'/V(v-v)- r+ (V-v)r-v = 
V-(r-vv) — %/2V((v-v)-r) — v-V; rev + 1/2(V -r)v-v 
(23) 


SadS\V XverXvt+V-ver-v) + 3/2f [doXverxXv+ 
do-vr-v] = '/2f'dSv-v (24) 

since 

v-do X (v X r) 

(25) 


doXvwvrxv=vXdovxXr= 


= v-rdo-v — v:vdo-r 


* —1/,.do-rv-v = '/2do Xv-r Xv — '/. do-vr-v 
(26) 
Energy Theorem 
The kinetic energy of a body of fluid is equal to the sum 
of the space integral of the dot product of the vector moment 
of the velocity and the rotation, of the space integral of the 
product of the scalar moment of the velocity and the diver- 
gence, and of one-half the corresponding surface expres- 
sions. 


The factor '/: in the surface terms is the consequence 
of the quadratic order of the terms. When a boundary 
is backed up by resting fluid and the integration ex- 
tended beyond it, the surface terms becoming space 
terms, the average velocity at the points of the former 
boundry is only one-half of the velocity at one side of 
the boundary, being the arithmetic means between that 
velocity and zero. Hence the same transition as before 
requires now an additional factor '/2, and the same is 
automatically obtained by the application of vector 
analysis. 

A theorem partially equivalent to the last one is dis- 
cussed in reference 2, Section 153. 


THE KuUTTA-JOUKOWSKI THEOREM 


The so-called Kutta-Joukowski theorem about the 
lift of a wing system moving with constant velocity in 
fluid otherwise at rest is another integral theorem, in 
that it gives only the total lift but not the distribution 
of the lift. The integrand differs from the terms pre- 
scribed by the basic vortex theorem in that the vorticity 
is multiplied by the constant velocity of advance, but 
not by the local flow velocity. 

This theorem can be generalized to apply to more 
general steady fluid motions, having vorticity. The 
entire space is divided into a specified region, for which 
the total force is to be computed, and into the remain- 
ing exterior space. Fictitious velocities within the speci- 
fied region are now computed by means of Biot Sav- 
art’s rule from the vortex elements of the exterior re- 
gion only. These velocities are substituted for the 
actual local velocities in the expression for the forces 
as specified by the basic vortex theorem. All these 
forces are then added up. The error committed by the 
use of the fictitious forces instead of the true local forces 
cancels out in the sum or integral, as will readily appear 
from the examination of Biot Savart’s rule, considering 
that in the force integral all pairs of vortex elements 
appear twice as products with opposite sign. The cor- 
rect total of the moments will be similarly obtained. 
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stant velocity of advance. 


Generalized Kutta Theorem 


The resultant of all forces necessary to act on a perfect 
fluid within a specified region to keep the motion steady 
can be computed by combining the velocities within that 
region, as determined from the exterior portion of the vor- 
tex system only, with the vortex system within the region, 
using otherwise the basic vortex theorem. 


VISCOSITY FORCES 


Concluding the discussion of some vortex theorems, 
it may be pointed out that the force per unit volume of 
the incompressible fluid ascribable to viscosity can also 


be simply expressed by means of the vorticity. This 
force is stated in most textbooks to be 
f=yV-Vv (27) 


that is, it is equal to the Laplacian of the velocity, 
multiplied by the viscosity modulus. See reference 2, 


AERONAUTICAL 


In Kutta’s case, these fictitious velocities are the con- 
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page 547, Eq. (4); reference 1, Vol. 3, page 49, Eq. 
(6.5). 

To transform this expression (27) into one containing 
the vorticity, there is combined the equation 


V:-Vv—-VV-v=V X(V Xv) (28) 
with 
V:-v=0 (28a) 
The viscosity force results then as 


re“2y x7 Xv =e2V Xe (29) 


Viscosity Force Theorem 


The resultant force per unit volume ascribable to the vis- 
cosity is equal to the rotation of the vorticity, multiplied by 
the viscosity modulus of the non-compressible fluid. 


It follows directly that this resultant force must 
vanish within the fluid in case of a potential flow. 

The modification of the Steady Flow Criterion by 
the inclusion of this viscosity term is easily suggested. 
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SUMMARY 


The analysis of a slender column of varying section is usually 
difficult. This paper suggests that the problem of the varying 
column be approached from the point of view of determining the 
variation of section which will give a certain previously desig- 
nated property. Three such properties are mentioned: 


(1) Constant M/EI along the length, under the critical load. 
(2) Constant Mc/EI under the critical load. 
(3) Minimum use of material. 

Cases (1) and (2) are represented by examples. Cases (1) and 


(2) come near to achieving (3). Case (1) is quite easy to apply. 





INTRODUCTION 


few TERM COLUMN as used herein refers to a slender 
compression member in which failure will occur by 
buckling. In the analysis of such a member, it is sup- 
posed that it is subjected to its critical load, that is, 
the load which will cause it to bend, but which it can 
support when in the bent condition. Under any smaller 
load, the member will remain straight. 

The logical reason for using a column of varying sec- 
tion, instead of a prismatic column, is to obtain certain 
properties not possessed by the latter. This being so, 
it appears that the most efficient method of dealing 
with a column of varying section is to decide what 
properties are desirable, and then to determine the 
variation of section which gives them; in many cases, 
this presents an easier problem than the analysis of a 
column of arbitrarily varying section. Properties 
which might be considered desirable in a column are: 

(1) Constant M/EI along the length, under the 
critical load. M, E and J represent bending moment, 
modulus of elasticity and moment of inertia of cross- 
section, respectively. This is an easy criterion to ap- 
ply. Also a buckling failure, for the type of member 
under consideration, is a bending failure. If M/EI is 
constant, the bending, as represented by the curvature, 
is the same at all points. 

(2) Constant Mc/EI under the critical load. Here 
c is the distance from the neutral axis of the cross-sec- 
tion to its extreme fiber. If this condition is met, the 
maximum stress will be the same at all cross-sections. 

(3) Minimum use of material. The first two condi- 
tions present comparatively easy problems. The third 
requires the use of Calculus of Variations. Blasiust 
has solved this problem for a few simple cases. In 
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some of them, the results are the same as those obtained 
by trying for properties (1) or (2). 


CENTRALLY LOADED COLUMN OF CONSTANT M/EI 
UNDER CRITICAL Loap P. No ENpD RESTRAINT 


In Fig. 1, consider first the centrally loaded column 
which is simply supported at the ends. For this case, 
M = —Py. It is customary in the analysis of slender 
columns to assume that the bending is similar to that 
of a beam, to the extent that the same differential equa- 
tion may be used, namely, d*y/dx*? = M/EI. Then 


d*y/dx* = —Py/EI 
Putting 
(1) 


a constant, d*y/dx? = —A, which integrates to y = 
(—Ax?/2) + Bx + C. Band C are constants of inte- 


Py/EI=A 


gration. From the condition that y = O when x = 
O,C =O. Using y = Owhenx = L 

O = (—AL?/2) + BL; B = AL/2 
Then 


y = (A/2) (Lx — x?) 


Eq. (1) can be rewritten as, y = AEI/P. Equating the 
two values for y, AEJ/P = (A/2) (Lx — x’). 


I = (P/2E)(Lx — x?) (2) 


This gives J as a function of x. The maximum J is re- 
quired at the center and is, J,, = PL*/8E. No restric- 
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tions on the shape of the cross-section have been im- 
posed in the above. The only requirement is that the 
bending axes of all sections be parallel. 

Blasius investigated the problem of the column of rec- 
tangular section of constant depth and varying width. 
He found that Eq. (2) gives the column of minimum 
volume for such a case, on the assumption that the 
bending is in one plane parallel to the depth. This as- 
sumption would not appear to be well founded, since 
the width is zero at the ends. 


CENTRALLY LOADED COLUMN OF CONSTANT Mc/EI 
UNDER CRITICAL LoaD P. No END RESTRAINT 


In this case, it is necessary to take the shape of cross- 
section into account. The case considered here is that 
in which the section remains of the same shape, circular 
or square, for instance, while changing in size. With 
this limitation, J = mc*, where m is a constant which 
depends on the shape; for a circular section, m = 7/4; 
for a square section, m = 4/3. From Fig. 1, 

d*y/dx* = M/EI = —Py/EI 
Thus, Mc/EI = —Pyc/EI. Putting Pyc/EI = A,a 
constant, and using J = mc*, Py/Emc® = A 


c = (Py/EmA)”8 (3) 
Also 

Py/EI = A/c 
and substituting the value of c from Eq. (3) 


Py/EI = (A‘mE/Py)”* = K*/y¥3 (4) 


where 
2 = (A‘mE/P)”3, a constant. 


Then 
d*y/dx? = —Py/EI = —K?/y"”* 


This can be integrated by use of the integrating factor, 
2(dy/dx). 
Thus, 


2(dy/dx) (d*y/dx*) = — 2(K*/y"*)(dy/dx) 
After integration 
(dy/dx)? = 3K*(d? — y”?) (5) 


in which the constant of integration appears as 3K°d?. 
A change of variables, letting y = 2°, leaves Eq. (5) as, 


otdeV d? — 2? = (K/V3)dx 
This integrates to 
— (¢/2)V d? — 2? + (d2/2) sin-! (¢/d) = (Kx/V3) + G 
Changing back to y, 


—(y¥3/2)V G2 — ys + (d?/2) sin (y'/3/d) = 
(Kx/V3) + G 


Since y = Owhen x = O,C, = O. Then, 
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x= sin~! is (6) 


3d2y/s 3 ys 
-¥i = ats 
4 K? 4 K? 2K 
Another boundary condition, y = O when x = L, gives 
d?x/2 = KL/V3; d2/K = 2L/V 3x 
From Eq. (4), 
y¥3/K? = EI/P 


Eq. (6) becomes 


L 3a? EI\”2 3x? ED 
i=-— = — —— —_ — —— 
® 4 =a 4 PL? 


The maximum [I is required when 


fs 3r? EI \Y4 7 3n? EI\"4 
sin -_ =) _ >? a ) = 1 





« 


So, 
Im = (4/3x?)(PL?/E) 


In terms of I,,, Eq. (7) is 


L vase oF I\“4 
es =F << — -__— in-!{— (8 
¥ al (5) _ + sin =) (8) 


Eq. (8) is sufficient for the determination of the cross- 
section at any point. By using Calculus of Variations, 
Blasius found that this result gives the circular column 
of minimum volume for a given critical load and length. 

It is interesting to compare this with the prismatic 
column. Consider first the maximum moment of in- 
ertia required for the same length and critical load. 
The coefficient of PL?/E is given below: 


0.1013 
0.125 
0.1352 


Prismatic column 
M/EI constant 
Mc/EI constant 


A saving in material results from using a column of 
varying section. The ratio of the volume of the varying 
column to that of the prismatic column for the same 


length and critical load appears below. Both figures 
are based on a circular section. 

M/EI constant 0.871 

Mc/EI constant 0.866 


From these figures, it appears that it makes little differ- 
ence which criterion is used, constant M/EI or con- 
stant Mc/EJ. The ease of application of the former is 
certainly a point in its favor. 


CENTRALLY LOADED COLUMN oF Constant M/EI 
UNDER CRITICAL Loap P. ENpDs RESTRAINED 


The entire column appears in Fig. 2a. The points of 
contraflexure, or points of zero moment, are marked 








COLUMNS OF VARYING 

















F, with the origin placed at the lower one. It is pos- 
sible to keep M/EJ, which equals d*y/dx*, of the same 
numerical value throughout the entire length, with a 
change of algebraic sign at each point F. For purposes 
of analysis, it is convenient to consider the column as 
divided into the three parts shown in Figs. 2) and 2c. 
The central part corresponds to the column discussed 
previously, with a change inlengthfrom LtokL. Then 
Eq. (2) may be applied and the moment of inertia is, 


I = (P/2E)(kLx — x?) (9) 


The maximum moment of inertia is J,, = Pk*L?/8E. 
While Eq. (9) gives true values only for the central part 
of the column, each end part, of length (Z/2)(1 — ), 
may be treated as a reflection about F, of an equal 
length of the central part. Thus, the required moment 
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of inertia at any point in the column may be obtained in 
terms of k. 

With a stable restraint at the end, My = n6, where @ 
is the amount of rotation at the end, and is a constant 
of proportionality. The use of makes it possible to 
express the end restraint quantitatively. For small 
deflections, 6 may be replaced by the slope, dy/dx. 
The general expression for the slope is, 


dy/dx = (A/2) (RL — 2x) 
At the lower end, 

dy/dx = (AL/2)(2k — 1) 
Thus 

Mo = (nAL/2)(2k — 1) 


Also, from Fig. 2c, Mo is equal to the product of P and 
the distance yo. In general, 


tila (A /2)(kLx — x?) 
At the lower end, 


— (A/2){(L2/4) — RL? + (3k°L2/4)] 


aa 
Then, 


My = —(PAL?/8)(1 — 4k + 3k?) 


Equating the two values for Mo, 


(nAL/2)(2k — 1) = —(PAL?/8)(1 — 4k + 3k?) 


n/PL = (1 — 3k)(1 — k)/(1 — 2k) (10) 


Eq. (10) is sufficient for the determination of k for a 
given 7. 

The moment of inertia equations developed above 
indicate zero dimensions for the cross-section at certain 
points. Slight modifications would of course be neces- 
sary at these points. 











The Measurement of Visibility 


MATTHEW LUCKIESH ano FRANK K. MOSS 
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ABSTRACT 


Certain fundamental principles involved in the measurement of 
visibility are discussed in relation to the use of the Luckiesh-Moss 
Visibility Meter. The construction, calibration and operation 
of this instrument are described. In view of the variety of 
visual problems in the aeronautical sciences and practices the 
present discussion is directed to general rather than specific tech- 
niques for investigating individual problems. 


VISIBILITY 


I AN OBJECT is visible, it is necessary that the size of 
the object, its contrast with its background, the 
brightness of the object or background, and the dura- 
tion of the image upon the retina each exceed certain 
individual threshold values.' * The threshold value of 
each of these variables is determined by the magnitudes 
of the other variables of the threshold stimulus. This 
is equivalent to stating that, within certain limits, a 
deficiency in one factor may be compensated by an 
augmentation of one or more of the other factors. It 
will be noted that visibility, as a measurable quantity, 
depends not only upon the physical characteristics of 
the object of regard but also upon the efficacy of the 
visual processes of the observer. However, if it be 
specified that all measurements of visibility are to be 
made by observers possessing normal vision or reduced 
to such a basis, the measured quantity then denotes a 
normal attribute of the object of regard expressed upon 
a rational scale. 

A common characteristic of present methods of meas- 
uring visibility is that such appraisals are based upon 
observations made under threshold conditions. This 
inherent restriction is due to the fact that the only 
unique visual stimulus between invisibility and maximal 
visibility is that of threshold visibility. Upon this 
basis, it is axiomatically assumed (1) that dissimilar 
objects of threshold visibility are equal in visibility and 
(2) that the same augmentation of one or more of the 
fundamental visual variables will cause these objects 
to be likewise equal in supra-threshold visibility. It 
is further assumed that the degree of supra-threshold 
visibility of a given object is proportional to the de- 
crease in the size, brightness, contrast or viewing-time 
required to reduce the object to threshold visibility.* 
For example, if the factors of brightness, contrast and 
time remain constant, the degree of supra-threshold 
visibility of a given object may be measured by the 
difference between its actual size and that size which 
would barely afford threshold visibility. This concept 
of the degree of supra-threshold visibility is analogous 
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The Luckiesh-Moss Visibility Meter bearing 


Fic. 1. 
scales calibrated in ‘‘Relative Visibility’? and ‘‘Relative 


Footcandles,’’ respectively. The latter scale is not dis- 
cussed in the present paper. Obviously various other 
scales of factors which influence visibility can replace the 
footcandle scale on the left. In fact, data for several 
such scales are already available. 


to the concept of ‘‘factor of safety,’’ as used in engineer- 
ing. This concept and a means of measuring supra- 
threshold visibility have been given little attention in 
a world in which things to be visible with certainty 
and with ease must be far above the threshold of visi- 


bility.” 
THE LucCKIESH-Moss VISIBILITY METER 


It follows from these basic principles that a ‘‘visibility 
meter’’ essentially involves (1) means for reducing the 
object of regard to threshold visibility and (2) means 
for expressing this reduction by impairment in the 
visual stimulus upon an appropriate scale. Thus a 
visibility meter may be calibrated in terms of a par- 
ticular fundamental visual variable while the requisite 
threshold conditions may be attained by altering any 
one or any combination of these variables. In the 
case of the meter developed by the authors, threshold 
conditions are obtained through the simultaneous re- 
duction of the factors of brightness and contrast. The 
practical advantage of this procedure is the enormous 
changes in visibility which are readily accomplished. 
However, the scale of visibility is calibrated in terms 
of the factor of size. The merit and usefulness of such 
a scale are largely due to the obviousness of the factor 
of size in seeing. In fact, size is more readily visualized 
than the other fundamental factors. 

The L-M Visibility Meter, illustrated in Fig. 1, con- 
sists essentially of two identical colorless circular gradi- 
ent filters, precise in optical density and diffusion, 
which may be rotated synchronously in front of the 
eyes of the observer while binocularly viewing an object 
while performing a visual task. One of these filters 
is illustrated in Fig.2. These filters incidentally reduce 
the apparent brightness of the visual field by absorption 
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One of the two gradient filters of the L-M 
The scale shown is calibrated in 
The 


Fic. 2. 
Visibility Meter. 
terms of threshold size as a measure of visibility. 
scale-values indicate size in minutes visual angle. 


but, of more importance, they reduce the degree of 
contrast between the object of regard and its back- 
ground, due to the slightly diffusing character (cf. “‘Q”’ 
factor in photography) of the photographic filters. 
Thus threshold conditions are established by rotating 
the gradient filters until the object of regard just 
emerges from, or disappears into invisibility; that is, 
As yet 


into the obscurity of the blurred visual field. 
it has been impracticable to produce filters which vary 
only the contrast and not the brightness but some 
promising possibilities are being investigated. 

In developing a quantitative scale of visibility, it is 
recognized that the specification of a unit of visibility 


is an arbitrary matter. However, this is equally true 
of such other units as the meter, gram, decibel, lambert, 
etc. The fact that a unit of visibility necessarily in- 
volves many arbitrary specifications does not lessen 
its practical usefulness. For example, the unit of 
visual acuity, as used in ophthalmology for rating 
human vision, involves all of the variables present in 
the specification of unit visibility as applied to objects. 
However, many of these variables disappear when 
visibility is measured upon a relative basis. This is 
likewise true in ocular refraction when the purpose of 
the observations is to measure refractive errors rather 
than to determine the ‘‘absolute’’ visual ability of the 
patient. 

The visibi‘ity meter is empirically calibrated by 
observing a series of test-objects of identical configura- 
tion but varying in size from the smallest visible to the 
largest encompassable by the human fovea.‘ The 
latter represents a practical upper limit in visibility 
since the effects of increasing beyond this value the size 
of a black object on a white background illuminated to 
a brightness of 10 footlamberts are practically negli- 
Such a series of test-objects is illustrated, with 


gible. 
In this case, the ‘“‘size’’ of 


some omissions, in Fig. 3. 


OF VISIBILITY 


< 3 a 5 6 


7 a 5 


10 20 


Fic. 3. Reproductions of the parallel-bar 
test-objects used in calibrating the L-M Visi- 
bility Meter. If these objects are viewed 
from a distance of 5 feet, the numerical 
values below the test-objects correspond to 
the visual angles in minutes subtended by the 
width of the bars or of the width of the space 
between the bars. 


the object is defined as the width of one of the parallel 
bars or the width of the space between them measured 
in minutes visual angle. It will be noted that the 
image of a “20-minute object’’ occupies a square 60 
minutes on each side, or practically the entire area of 
the fovea. These test-objects were repeatedly viewed 
by a group of observers possessing normal vision and 
the mean position of the gradient filter at which each 
test-object became visible was marked upon the visibil- 
ity scale (Fig. 2) attached to the filter. 

In principle, the threshold visibility of any object, 
viewed under any condition by any observer is equiva- 
lent to one of the infinite series of standard test-objects. 
A value of “‘1”’ on the scale of visibility indicates that a 
standard test-object, whose critical details subtend a 
visual angle of one minute, is of threshold visibility 
when viewed through a region of the filters of the meter 
so designated. Hence any object which is barely 
visible through the same portion of the filters is assigned 
a value of unity in visibility. Specifically, the visibility 
of such an object, under the conditions of measurement, 
is equivalent to that of a 1-minute standard test-object 
as viewed by an observer possessing normal vision 
and under a brightness-level of 10 footlamberts. The 
latter object is practically equivalent in visibility to 
the ‘‘20/20” characters (assumed standard for normal 
vision) on the test-chart of the ophthalmologist. 
Similarly, a scale value of ‘‘2’’ indicates that the visibil- 
ity of the object of regard is equal to that of a 2-minute 
standard test-object or to that of the ‘‘20/40” char- 
acters of the test-chart, etc. Obviously, these values 
of visibility are significant in an absolute sense, and 
they are uniquely related to the universally assumed 
standard of normal vision, only when the observations 
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are made under the same conditions as those prevailing 
during the calibration of the meter. 

Since a scale value of ‘‘1”’ represents threshold visibil- 
ity under the specified conditions of calibration, it is 
obvious that higher values represent supra-threshold 

For example, if the visibility of a given 
object is ‘‘5,’’ such an object is equivalent in visibility 
to a standard parallel-bar object five times as large as 
the smallest visible under the same conditions of ob- 
servation. Thus the object of regard may be said to 
possess a visual ‘factor of safety”’ of five with respect to 
the variable of size of an equivalent standard object. 
However, in the interpretation of such data it should 
be recognized that the factor of size obeys a “law of 
diminishing returns” as the object becomes larger. In 
consideration of this fact another scale of visibility has 
been developed? in which the units of visibility through- 
out its entire range are equivalent in visual effective- 


visibilities. 


ness. 


VISIBILITY MEASUREMENTS 


The visibility meter may be used in practice in two 
different ways: (1) as a means for measuring visibility 
in terms of a rational absolute scale; or (2) as a means 
for appraising the visibilities of two or more objects 
upon a relative basis. If significance is to be attached 
to the scale values in an absolute sense it is obvious that 
all measurements must be made under the same physical 
conditions as those prevailing during the calibration of 
the meter. These conditions may be reproduced 
readily. It is also necessary that the observers shall 
possess normal vision and the same introspective con- 
ception of threshold visibility as the observers used in 
the calibration of the meter. Since these human 
variables are usually beyond control, it becomes neces- 
sary to obtain ‘correction factors’ for the individual 
observers. For example, if a 10-minute standard object 
is found by an individual observer, to possess a visi- 
bility rating of “‘8,’’ all measurements taken by this 
observer should be increased by 25 percent. Obviously, 
such corrections would be required in meteorological 
surveys conducted by different observers. 

The simplest and most general use of the visibility 
meter is in determining the relative visibilities of various 
objects. When used for this purpose, the only require- 
ment is that the same conditions prevail throughout a 
given series of measurements. Thus the measurements 
may be made at any distance, under any level of il- 
lumination or lighting condition and by any observer. 
Furthermore, the values of relative visibility so ob- 
tained may be reproduced with considerable accuracy 
by another observer and under other conditions for 


observation. For example, a group of children in the 


sixth grade of a Sight-Saving class obtained practically 
the same ratio (within less than 2 per cent) in visibility 
between two typographical samples as did a group 
of technically trained adult subjects possessing normal 
or near-normal vision. 
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It should be recognized that a single measurement 
whose magnitude may be markedly affected by human 
vagaries possesses but little intrinsic reliability. Hence 
only average values, obtained from an adequate number 
of measurements in each case, should be given quantita- 
tive consideration. Obviously, the number of observa- 
tions required to produce a reliable measure of visibility 
depends upon the skill of the observer and upon the 
significance to be attached to the value found. To 
discriminate between objects which differ but little in 
visibility, it is necessary to make more observations 
than are required where the differences are large. In 
general, it is advisable to test by statistical methods the 
reliability of all measurements of visibility. For ex- 
ample, if the difference in visibility between two objects 
is found to be several times as large as its probable 
error, the difference may be considered to be significant. 
As experience is gained in the use of the meter, and in 
interpreting visibility data, there will be less need for 
detailed statistical analyses of the data. Analyses of 
the precision obtainable with the L-M Visibility Meter, 
including a study of the variability within a single 
series of measurements, as well as the variability among 
average values obtained on different occasions, are 
available elsewhere.* 


APPLICATIONS IN AERONAUTICAL SCIENCES 


The visibility meter should be of value in many ways 
in aeronautics such as (1) in appraising the influence of 
fog, smoke, etc., upon the visibility of appropriately 
placed standardized targets, (2) in measuring the loss 
in visibility due to transparent media, windows, goggles, 
etc., in the line of vision, (3) in determining the visibil- 
ity of various signals, different types of instrument dials, 
etc., and (4) in appraising the visual abilities of pilots 
and observers to distinguish low contrasts at night as 
well as during the day. Some of these purposes require 
absolute measurements involving the observation of 
standardized test-objects under rigidly specified condi- 
tions, while others require merely relative measure- 
ments under any constant set of conditions. An ap- 
preciation of the principles involved in the measure- 
ment of visibility should provide an adequate basis re- 
vealing uses for the meter and for developing tech- 
niques of measurement in specific fields of application. 
Obviously knowledge in the fields of seeing and of aids 
to seeing can only be established by meaningful meas- 
urements as in other fields. 
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The Transmission of Light in the Atmosphere 
with Applications to Aviation 
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ABSTRACT 


The transmission of light through the atmosphere is determined 
by the absorption and scattering which it suffers. It is shown 
that absorption is negligible in the visible spectrum but becomes 
increasingly strong in the infrared. Scattering is due to sus- 
pended particles ranging in size from the air molecules to rain 
drops. The scattering coefficient changes with the wave length, 
or color, of the light only when the scattering particles are smaller 
than the wave length. Since the suspensoids which prodtce 
low visibility are large compared to the wave length of visible 
light and the near infrared, no advantage results from the use of 
colored lights or the infrared in such cases. Photography in the 
near infrared is useful in cases when the visibility is limited by 
small particle haze. 

Theory shows that the visual range is determined largely by 
the scattering coefficient of the atmosphere. Since it is often 
difficult to select suitable marks for the estimation of the visual 
range, it is suggested that an instrument be developed for deter- 
mining the scattering coefficient. 


INTRODUCTION 


Fyre THE transmission of light through the 
atmosphere has long been of interest to astrono- 
mers and other workers in the natural sciences, the ad- 
vent of air transportation has served to emphasize the 


practical importance of this subject. In spite of the 
great advances in navigation instruments, pilots still 
depend to a considerable extent on visual observation 
of prominent landmarks by day and of lights and 
beacons at night. The distance at which such marks 
can be seen is largely determined by the transparency 
of the atmosphere. The same factor is also of major 
importance for aerial photography and for military 
reconnaissance and bombing missions. In an attempt 
to extend the useful range of observation, lights of vari- 
ous colors and instruments utilizing infrared radiation 
have been tried with varying degrees of success. 

It is the purpose of this paper to discuss the prin- 
ciples governing the transmission of light in the atmos- 
phere and to apply them to some of the problems of 
special interest for aviation. In particular the factors 
governing the visibility of objects by day and lights at 
night will be discussed and also the merits of infrared 
or other invisible radiations for signalling and photog- 
raphy. No attempt will be made to present a com- 
plete discussion of the theory of the visual range nor 
of the details of photography and special instruments, 
since this is available elsewhere. 
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THE DIMINUTION OF LIGHT IN THE ATMOSPHERE 


The most useful measure of light for the purposes of 
this paper js the flux density which is flux per unit area. 
Flux is a measure of the visible power, but the concept 
of flux density may also be used in the case of invisible 
radiation, in which case the flux is considered to be the 
total radiant power. 

If a beam of light of restricted angular extent is di- 
rected through the atmosphere, the flux density of the 
beam will decrease with distance from the source due to 
three effects. The first of these is the inverse square 
law which states that the flux density will decrease 
inversely as the square of the distance from the source. 
This decrease is due solely to the geometry of the optical 
system. Second, the flux density will decrease due 
to the absorption of the light by the atmospheric gases 
and suspended matter. Finally, there will be a decrease 
due to scattering by the air molecules and the various 
atmospheric suspensoids. In this paper the term scat- 
tering is taken to mean any process which changes the 
direction or spatial distribution of the light and hence 
includes reflection, refraction and diffraction. Al- 
though scattering does not involve any energy trans- 
formations, it is the principal factor governing the trans- 
mission of visible light in the atmosphere. 


ABSORPTION IN THE ATMOSPHERE 


The absorption due to the fixed gases of the atmos- 
phere is entirely negligible in the visible and infrared 
regions of interest to this discussion. Oxygen in both 
the atomic and molecular forms and ozone absorb 
strqngly in the ultraviolet and combine to set a lower 
limit of the solar spectrum at about 0.34, but this is not 
of much importance in the present investigation. 
The important absorbing gases in the atmosphere are 
carbon dioxide and water vapor. Carbon dioxide has 
three principal absorption bands, one at about 2.7y, 
a second at 4.5u and the third and most intense ex- 
tending from about 12 to 16. Carbon dioxide ab- 
sorption thus becomes of significance only in the infra- 
ted. Water vapor is also quite transparent within the 
visible spectrum (.4 to .74) but absorbs in several 
regions in the near infrared and has a complex and in- 
tense absorption spectrum in the region beyond 5x. 
The general features of this portion of the spectrum are 
intense absorption from about 5 to 6y, a relatively trans- 
parent region from about 9 to 15u and rapidly increas- 
ing absorption for wave lengths greater than 20. 
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Certain suspensoids, such as smoke and dust, may 
absorb visible light to some extent, but the scattering 
caused by such particles appears to be large compared 
to the absorption. Water drops which occur in the at- 
mosphere as fog, clouds, drizzle, rain, etc., are essen- 
tially transparent in the visible spectrum. However, 
the absorption of water increases rapidly in the infrared, 
a layer 1 mm. thick being opaque to wave lengths 
greater than 2.4u. Because of the effects of scattering, 
it cannot be assumed that a cloud absorbs the same 
fraction of the radiation as an equivalent sheet of water, 
but experiments indicate that an average cloud 100 
meters or more in thickness acts essentially,as a black 
absorber in the infrared. 


SCATTERING IN THE ATMOSPHERE 


It is clear from the above discussion that absorption 
in the atmosphere is important only in the infrared re- 
gion. The major factor which determines the trans- 
mission of visible light in the atmosphere is scattering 
by the various suspended particles. In approximate 
order of increasing size of the particles the important 
scattering elements in the atmosphere are the air mole- 
cules themselves, smoke, haze, dust, cloud and fog par- 
ticles and the various precipitation elements. 

In the case of a parallel beam, the reduction in flux 
density due to scattering as a function of distance is 
given by the well-known exponential law: 


E/E, = e~™ (1) 


where £ is the flux density at the distance x along the 
beam 
E, is the flux density for x = 0 
k is the scattering coefficient 


Since parallel beams are not obtainable in practice, the 
inverse-square law must be combined with the above 
exponential expression, giving 
IR —kx j9 ‘< 
E/E, = Toe So (2) 


where J) is the candlepower of the source, or the beam 
candlepower in the case of a collimated source. 

In general, the scattering coefficient is a complex 
function of the ratio of the particle diameter and the 
wave length of the light. For the special case of par- 
ticles small compared to the wave length Rayleigh’ 
showed that the scattering coefficient varied inversely as 
the fourth power of the wave length. This relation 
has often been used erroneously to predict the scatter- 
ing due to particles which are not small compared to 
the wave length. Table 1 gives the approximate 
range of ‘“‘diameters’’ of the several types of scattering 
particles in microns. The visible spectrum extends 
from about .4 to .7u and the near infrared from 17 to 
20u, according to Strong’s? classification. It is clear 
from the table that for visible light, Rayleigh’s in- 
verse fourth power relation is valid only for scattering 
due to molecules and the smaller haze, dust and smoke 
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TABLE | 





Type of Particle Approximate Diameter of Particles 


molecules 10-4 microns 
smoke .01-.5 

haze .1-10 

dust . 1-50 
cloud 5-50 

fog 5-100 
drizzle 107-5 X 10? 
rain 5 & 10*-7 X 10° 


particles. It is also evident that it would be necessary 
to go far into the infrared to obtain a wave length large 
compared to fog and cloud particles. In this region the 
absorption of water vapor becomes very intense so that 
the reduction in scattering would probably not result 
in an increased transmission. 

The variation of the scattering coefficient when the 
particles are not small compared to the wave length 
has been treated by Stratton and Houghton® for the 
case of spherical water drops in air. The scattering 
coefficient was found to be rur*K where r is the radius 
of the drops, m the number of drops per unit volume of 
air and K is a function of 27r/X as shown in Fig. 1. 
In the original paper the scattering coefficient was 
given as 2rnur*K, but it now appears, as the result of 
some experimental work, that the factor 2 should not 
be included. This disagreement is not due to a numeri- 
cal error but to something more fundamental which is 
not yet completely understood. 








20 — : ——, 
| | 























| | 
| 
| j 
| Ape ERS 
. H 
1@) 10 20 30 40 50 
2 
m% 
Fic. 1. Curve of K, the ratio of the scattering coeffi- 


cient to the projected area of the drops for spherical water 
drops in air, versus 2nr/x. 


The portion of the curve of Fig. 1 for 27r/\ less than 
about unity corresponds to the region of Rayleigh 
scattering. For large values of 27r/A, exceeding about 
30, K is approximately unity. It is clear that a rapid 
decrease in scattering with increasing wave lengths can 
be expected only along the steep portion of the curve 
to the left of the first maximum, that is for wave 
lengths exceeding the radius of the particles. With the 
aid of Table 1 the minimum wave length for each type 
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of suspensoid may be determined. It is clear that no 
great reduction in the scattering produced by particles 
as large as fog and cloud drops can be expected in the 
visible spectrum or in the near infrared. This result 
has been confirmed by recent measurements of Sander- 
son.* 

Although Fig. 1 applies strictly only to spherical, non- 
absorbing water drops, it is believed that it is also ap- 
proximately correct, at least in the regions to the left 
of the first maximum and beyond the second maximum, 
for non-spherical particles of approximately the same 
index of refraction. It will be noted that as K ap- 
proaches unity the scattering coefficient rur*K ap- 
proaches the total projected area of the drops in a unit 
volume of air. That is, all of the light intercepted by 
the drops is scattered. It would appear from geo- 
metrical optics that this would hold for non-absorbing 
particles of any shape and index of refraction as long as 
the particles are large compared to the wave length, 
although the spatial distribution of the scattered light 
will be different in each case. If the particles also 
absorb, it is not proper to add the scattering and absorp- 
tion coefficients. If the absorbing particles are large 
compared to the wave length, it would appear that the 
total of the absorption plus scattering could not exceed 
the scattering due to a non-absorbing particle of the 
same size. It is not possible to predict the fractions 
scattered and absorbed except in special cases. The 
case of absorbing particles which are not large compared 
to the wave length is very complex and no complete 
solution to this problem has been found. 


APPLICATIONS TO LIGHT SIGNALS 


. 


In the case of visual signal lights or beacons, it is 
clear from the above general discussion that no change 
in transmission with wave length can be expected unless 
the scattering particles are small compared to the wave 
length. Low visibility is generally due to particles 
larger than the wave length of visible light, and experi- 
ence as well as theory has demonstrated that colored 
lights are of no advantage in such cases. Occasionally 
a fine-particle haze may reduce the visibility to from 
three to five miles. If it is assumed that Rayleigh’s in- 
verse fourth power law applies, a relatively small change 
in wave length will produce a considerable change in 
transmission. However, the sensitivity of the eye, 
which is a maximum at about .55 also varies rapidly 
with wave length. These two effects will combine 
to give a maximum visibility in the neighborhood of 
0.64, which corresponds approximately to light from a 
sodium lamp. However, this conclusion applies only to 
the special case of a haze composed of very small par- 
ticles and, in general, no appreciable gain can be ex- 
pected as has been shown experimentally by Luckiesh 
and Holladay.’ In particular, color filters over white 
lights definitely reduce the visibility of the light, since a 
filter decreases the intensity of the source even in the 


pass band of the filter. Colored or infrared lights may 
sometimes be necessary for purposes of identification or 
secrecy, but no appreciable gain in visibility is produced 
even when the lights are viewed against a general back- 
ground illumination as shown by Langmuir and Westen- 
dorp.® These authors also found that the visibility 
of flashing point sources never exceeds that of a steady 
light of the same intensity. 

As shown in Eq. (2), the flux density at a distance x 
from a point source of intensity J is directly propor- 
tional to the intensity and to the transmission of the 
atmosphere and inversely proportional to the square of 
the distance. The light will be visible if the flux den- 
sity at the eye exceeds a certain minimum value which 
may be taken as about 2 X 10~' lumens per sq. cm.’ 
The candlepower required to make the light visible at a 
distance x is then 


I = Eyx*e™* (3) 


where EF is the threshold flux density for the eye. As 
will be indicated below, the daytime visual range of @ 
dark object against the horizon sky is given by 


Xm = (log, 1/e)/k (4) 
where ¢ is the threshold of contrast which has a value 
of about .02. Thus, with the aid of Eq. (4), Eq. (3) 
may be written 

I = Equ?/(e)*/*™ (5) 


As an example of the application of this formula suppose 
that the transmission of the atmosphere is such that 
the daytime visual range would be one mile. It is de- 
sired to find the candlepower of a light required to make 
it visible at two miles. 


I [2 x 10-" X (3.22)? & 10?°] /(.02)? 
= 5190 c.p. 


II 


The only practicable means for increasing the visual 
range of a light is to increase its candlepower in ac- 
cordance with Eq. (5). As shown by the above ex- 
ample, it is relatively simple to provide lights which 
can be seen at night at distances considerable greater 
than the daytime visual range in the same atmosphere. 


APPLICATIONS TO PHOTOGRAPHY 


The availability of films which are sensitive in the 
near infrared has made it possible to take advantage 
of the fact that on days with relatively good visibility 
the scattering particles are usually small compared 
to the wave length of visible light. In such cases, the 
atmosphere may be several times as transparent in the 
near infrared as in the visible. It also follows that the 
veiling glare due to the scattered light will be corre- 
spondingly reduced in the infrared. These factors 
combine to yield much clearer photographs in the infra- 
red. Since infrared photography is limited to the very 
near infrared (limit about 1.2), it will not permit pho- 
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tographs to be taken through fog or other suspensoids 
in which the particles are not less than about 0.5y in 
diameter. Unless distant objects appear bluish, it is 
extremely doubtful if infrared photography will be ad- 
vantageous. For a detailed presentation of the many 
problems involved in infrared photography, the reader 
is referred to the recent book by Greenwood.* 


DAYTIME VISIBILITY 


In the absence of color contrast, an object is visible 
when the ratio of the difference in brightness between 
the object and its background to the background bright- 
ness exceeds a certain minimum value. The value of 
this brightness threshold is .01 to .02 under normal con- 
ditions, although it is subject to considerable variation. 
As the distance between an observer and the object 
increases, the apparent brightness of a dark object 
increases due to the sun and sky light which is scattered 
into the eyes of the observer by the atmospheric sus- 
pensoids in the line of sight. A distance will finally be 
reached at which the apparent brightness of the object 
will approach that of the background so closely that 
the eye cannot distinguish the object. This limiting 
distance which is called the ‘‘visual range’ or the “‘visi- 
bility” is a function of the scattering coefficient of the 
atmosphere and of the reflectivity of the object and the 
background. Expressions for the visual range for a 
number of cases have been derived but only the re- 
sults will be presented here. For more details the 
reader is referred to the excellent book by Middleton.’ 


Case 1. A black object viewed against the horizon sky. 
Xm = (log, 1/e)/k 
where x,, = the visual range 
€ = contrast threshold (.01 to .02) 
k = scattering coefficient 
Case 2. An object of albedo* A viewed against the 
horizon for a totally overcast sky. 
%m = (1/k) log. (1 — A/2)/e 
Case 3. Object of albedo A; viewed against a back- 


ground of albedo Az. Totally overcast sky. 
Object and background at the same distance. 
Xm = (1/k) log, (1/2e) [((A1 — Az) + 


2e — €A,] 


If (A; — Az) is reasonably large, the latter 
two terms in the brackets may be neglected. 


An inspection of the above formulas shows that in 
each case the visual range is inversely proportional to 
the scattering coefficient. The expression given above 
for Case 1 has been experimentally verified by Hul- 
burt.” 

The ratio of the visual range of an object of albedo A 
to that of a black object when both are seen against the 
horizon sky under the same conditions is 


* Albedo = diffuse reflection coefficient. 


log { [1 — (A/2)]/e}/log (1/e) 


The minimum value of this ratio is .82, which corre- 
sponds to the maximum possible albedo of unity. It 
follows that as long as objects are viewed against the 
horizon sky, the albedo of the object has little effect 
on the observed visual range, although dark objects 
should be selected if possible. When the object is 
viewed against a background of albedo A, (Case 3), 
the visual range is almost directly proportional to the 
log of the difference in albedos of the object and the 
background. It is clearly possible for the visual range 
of an object seen against a terrestrial background to be 
considerably less than the visual range of the same ob- 
ject seen against the horizon sky. The maximum pos- 
sible value of (A; — Ag) is unity, and for this case the 
visual range is again 82 per cent of the visual range of 
a black object against the horizon sky. This is of con- 
siderable importance in aviation, since visibility esti- 
mates are usually based on observations of objects seen 
against the sky while a pilot views objects against the 
ground. 


If the scattering coefficient k could be assumed con- 
stant in all directions, it would be possible to obtain a 
relation between the visual range of a black object seen 
against the sky and an object of known albedo viewed 
against the ground. The suspensoids which determine 
the scattering coefficient in the lower atmosphere are 
usually stratified horizontally and decrease rapidly in 
number with altitude. This is borne out by the ex- 
ceptional horizontal visual ranges often observed 
at high levels. Thus the mean scattering coefficient 
in the vertical direction is commonly considerably less 
than in the horizontal direction at the ground. This 
difference is often so marked that a well-defined ‘‘cone 
of visibility,”’ with its apex at the object, is formed when 
an object on the ground is viewed from the air. The in- 
cluded angle of the cone will usually decrease as the 
elevation of the point of observation is increased. 


Airplanes are often visible from the ground when the 
pilot is unable to observe any objects on the ground. 
As often suggested, this may be due to reflection of sun- 
light from the sharp “top” of a haze layer which may 
form at a temperature inversion. However, even in the 
absence of this phenomenon, the visual range of an 
airplane will often exceed that of terrestrial objects 
because of the difference in background (Case 2 vs. 
Case 3). This would always be true except that the 
contrast threshold increases as the angle subtended by 
the object at the eye decreases and prominent terrestrial 
objects are usually larger than airplanes. 


In connection with camouflage problems, it should be 
noted that an object can never appear as bright as the 
sky unless it is reflecting sunlight. Furthermore, a 
“black” airplane would be visible about 20 per cent 
further away than one having a diffusely reflecting sur- 
face of albedo unity. 
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MEASUREMENT OF THE VISUAL RANGE 


It is customary to estimate the visual range by noting 
the most distant prominent object which can be seen. 
It is seldom possible to select a series of marks located 
at appropriate distances which are all seen against the 
horizon sky, are all nearly black and of sufficient size. 
These difficulties could be overcome by the use of a 
suitable instrumental technique. This would have the 
further advantage of making measurements at different 
stations directly comparable. 

The most useful quantity to measure is the absorption 
coefficient, since it would then be possible to interpret 
this in terms of the visual range either of lights of 
known candlepower or of objects seen either against the 
sky or the ground. It would also be desirable to know 
the variation of the absorption coefficient with eleva- 
tion so that the oblique visual range could be estimated, 
but an accurate determination at the ground level 
should be attempted first. 

The scattering coefficient may be determined most 
directly by measuring the transmission of light over a 
path of known length,'! but the equipment is too expen- 
sive for general use. A number of empirical meters 
have been developed in Europe in which ground glass 
discs of various degrees of roughness are interposed be- 
tween the eye and the object until the object is no longer 
visible.'!* These instruments do not measure the scat- 
tering coefficient and they depend to some extent on the 
human equation. In spite of these limitations, they 
have given good results but have not been used in this 
country. It is suggested that a more suitable instru- 
ment would be a visual photometer adapted to measur- 
ing the difference in brightness of a relatively dark 
object at a known distance and the horizon sky. Such 
a measurement would determine the scattering coeffi- 
cient directly. Goldberg’ has invented a photometer 
which could be used for this purpose, and it is believed 
that a simpler and inexpensive photometer of this type 
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might be developed. Such an instrument could also be 
used at night if arranged so that a distant light of known 
candlepower could be compared with a standard light 
in the instrument as was done by Middleton.’ 
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An Advanced Method for the Determination 
of Intermediate Wing Contours 


T. GARDNER HILL* anp MAXWELL B. BASSETT? 
The Glenn L. Martin Company 


SUMMARY 


As speeds are approaching critical values and the use of twisted 
and dis-similar airfoil sections is becoming more frequent, the 
accuracy of the full scale master lines from which wings are 
constructed is becoming more critical. A geometrical method is 
suggested for the determination of intermediate wing contours for 
wings with dis-similar or twisted airfoils. 

The problem, briefly stated, is to determine intermediate 
wing contours when the root and tip airfoil sections and their rela- 
tion to each other are given. It is assumed that the most de- 
sirable wing aerodynamically and structurally, where there is no 
reverse of curvature in the limiting airfoils, is a straight line gen- 
erated surface about these two limiting sections. (Where there 
is reverse of curvature in just one of the airfoils, a compromise is 
desirable. ) 


LIMITATIONS AND ERRORS OF CUSTOMARY METHOD 


_ CUSTOMARY method of determining intermedi- 
ate contours involves the selection of points on the 
limiting airfoils at corresponding percentages of the 
airfoil chords and the drawing of straight lines between 
these points, which straight lines are then divided pro- 
portionally to the wing rib stations to give points for 
the drawing of the intermediate contours (see Fig. 1). 


r INTERMEDIATE CONTOUR POINTS 
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This method arose from the much used straight line 
tapered wing having similar airfoils with chords in a 
common plane. For this special case the percentage 
lines are straight lines on the generated surface. How- 
ever, as soon as the limiting airfoil sections are made dis- 
similar or are twisted with reference to each other, then 
the percentage lines are no longer the location of the 
generating lines and the percentage points on the gen- 
erated surface form curved lines. So far, the error 
involved in assuming the percentage lines as straight 
lines and forcing the wing skin to fit the contours thus 
determined has not proved of practical importance. 
Nevertheless, wing ribs made from contours determined 
by this method will force a series of flats or reverses in 
fairing in the wing skin and as the airfoil dis-similarity 
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and wing twist increase, this forcing will reach a struc- 
tural limit, especially with corrugated wing box con- 
struction. Furthermore, it is believed that even very 
slight surface irregularity will be aerodynamically pro- 
hibitive in the future. 


TANGENTIAL METHOD OF CONTOUR FAIRING 


A uniform sheet wrapped around any two limiting 
smooth contours will have straight lines on its surface 
in certain positions only. Fig. 2 illustrates this fact; 
the dotted lines on this figure connect the conventional 
percentage points, the full lines show the location of 
straight lines as determined by the use of a straight 
edge. In this case it is apparent that the assumption 
of straight lines along the percentage points is in error. 








To determine the location of the straight lines (gen- 
eratrices) on our surface let us place the limiting air- 
foils in parallel planes A and B and lay our surface on a 
third plane D (see Fig. 3). It becomes apparent that 
the surface makes line contact with plane D through 
points a and b on the two limiting airfoils where the air- 
foils are tangent to the plane. Since these tangents lie 
in a common plane in one projection and in parallel 
planes in another projection, then in any other projec- 
tion they will appear as parallel lines. Thus, the cri- 
terion for the location of the generatrices is that they 
connect points of parallel tangency on the limiting 
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curves. For a more rigorous proof of this conclusion 
see ref. 1. 

As an example, Fig. 4 is an exaggerated wing with 
dis-similar and twisted limiting airfoils A and B. By 
constructing parallel tangents such as ab and cd the 
generatrices are readily located around airfoil A to 
point /. This leaves area /nm for which we have points 
on airfoil B but no corresponding points on airfoil A. 





In the practical case of an actual skin wrapped around 
airfoils A and B the surface /mn would contain straight 
lines in some such locations as gw, hs, it, etc., and would 
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terminate in a curved trailing edge. However, for 
structural reasons a straight trailing edge is desirable 
and so it is best to consider the area /nm as a portion of 
a cone with apex at / and draw generatrices of the cone 
as shown by the dotted lines in Fig. 4. 

The envelope of the generatrices should fall on a 
smooth curve called the cuspidal (Q, Fig. 4). The cus- 
pidal curve is useful as a check on the proper location 
of the generatrices and there should be such a curve 
asymptotically approaching each generatrix where the 
rate of curvature of the two limiting curves becomes 
equal. The location of this asymptote may be of 
aerodynamic use in determining the location of points 
of equal pressure on the wing surface. 


CONCLUSIONS 


1. The percentage line method for the determination 
of intermediate wing contours is accurate only for 
straight taper wings with similar limiting airfoils having 
chords in a common plane. 

2. For the determination of correct intermediate 
contours on wings having dis-similar limiting airfoils or 
wing twist, the following procedure is suggested: 

(a) Lay out the two limiting airfoil sections to 
full scale and in their correct angular relation. 

(6) Through arbitrary points on the limiting airfoils 
draw generatrices by the tangential method. 

(c) Divide the generatrices proportionally to the 
spanwise locations of the intermediate contours desired 
and draw intermediate contours through these points. 

3. For the application of the tangential method to 
control surfaces, see ref. 2. 
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Sn WRIGHT BROTHERS Lecture, commemorating 
the thirty-eighth anniversary of the first flights 
made by Orville and Wilbur Wright at Kitty Hawk, 
was presented by the Institute on Wednesday, Decem- 
ber 17, 1941, at Columbia University, New York City. 
This was the fifth of these annual lectures which are 
endowed by the Vernon Lynch Fund established in 1937 
by the late Edmund C. Lynch in memory of his brother. 

The lecturer was Richard V. Southwell, Professor of 
Engineering Science at Oxford University, England. 
Professor Southwell holds an M.A. degree from Oxford 
University and also from Cambridge University where 
he was one time lecturer in Mathematics. He is a 
world renowned specialist in aerodynamics, hydro- 
dynamics, theory of structures, elasticity and mathe- 
matical applications to aeronautical problems. 


DINNER IN HONOR OF THE LECTURER 


Presentation of the lecture was preceded by a dinner 
in honor of Professor Southwell, held at 6:30 p.m. at the 
Men’s Faculty Club of Columbia University. It was 
attended by 200 Institute members and their guests. 

Frank W. Caldwell, President of the Institute, acted 
as Toastmaster. 

Mr. Caldwell introduced several of the distinguished 
guests. In presenting the guest of honor, he said: 

‘Professor Southwell, I wish to express our apprecia- 
tion of your making this trip to America to address the 
Institute on the occasion of the Wright Brothers Lec- 
ture. I also wish to welcome you to America. We had 
expected to welcome you as a distinguished British scien- 
tist, as a representative of British aeronautical science 
and of the Royal Aeronautical Society. Tonight, we 
can go further. As American citizens, we can welcome 
you as a representative of a staunch ally.”’ 


PRESENTATION OF THE LECTURE 


The lecture was held at 8:30 p.m. before a capacity 
audience in the main lecture theatre of the Pupin 
Physics Laboratories Building of Columbia University. 
Chairman of the meeting was Dr. J. C. Hunsaker, 
Chairman of the National Advisory Committee for 
Aeronautics and Past President of the Institute. 
Honorary Chairmen were C. R. Fairey, Deputy Direc- 
tor-General of the British Air Commission in Washing- 
ton and Past President of the Royal Aeronautical 
Society of Great Britain, and Dr. Charles Darwin, 
Chairman of the British Central Scientific Office in 
Washington and Director of the National Physical 
Laboratory, England. Mr. Caldwell opened the meet- 
ing and introduced Dr. Hunsaker. 

Dr. Hunsaker: “We meet tonight on the 38th 


anniversary of the Wright Brothers’ first flight over 
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Prof. Richard V. Southwell 


the beach at Kitty Hawk, North Carolina, in 1903, 
and for the fifth Wright Brothers Lecture. To mark 
this event, and the most important invention of the 
century, the Royal Aeronautical Society of Great 
Britain has held for many years at London the Wilbur 
Wright Memorial Lecture. Some twenty years ago, 
I was honored by an invitation to give this lecture, 
and was most hospitably and kindly helped through 
what seemed to me a pretty difficult assignment. 

“IT had hoped to welcome to this meeting as an 
Honorary Chairman, Mr. C. R. Fairey, Past President 
of the Royal Aeronautical Society and a member of 
the Council of the oldest and most distinguished aero- 
nautical society in the world. Mr. Fairey is in our 
country in connection with the joint effort to create 
fleets of airplanes adequate to our grave needs and was 
at the last moment prevented from coming. However, 
I have a cablegram from Mr. Griffith Brewer, who is so 
ably leading the activities of the Royal Aeronautical 
Society in London today: 


“On behalf of my council and all members, I convey 
to you on the occasion of the reading of the Wright 
Brothers Lecture by our distinguished member, Professor 
Southwell, our best wishes for a successful meeting. 
On this particular occasion I would add how profoundly 
we are all moved by the fact that members of the two 
societies are now in tt all the way until the freedom of 
mankind and peaceful flight are assured to us all. 

GRIFFITH BREWER, President 
The Royal Aeronautical Society 


“‘T now wish to call upon Lester D. Gardner, Execu- . 
tive Vice-President of the Institute, who has some an- 
nouncements to make.”’ 
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Mr. Gardner: “At this first meeting of Institute 
members since the passing on December 9 of our Bene- 
factor Member, Hart O. Berg, I want to take the 
opportunity to express publicly the great personal loss 
that is felt by the Officers and Council of the Institute 
and many members who knew him personally. He 
will be missed at our meetings where his gracious 
courtesy was always in evidence. Hart O. Berg was 
the one who really started the idea of the Institute 
having an Aeronautical Archives through his gift of a 
priceless collection of papers, pamphlets, books of 
clippings and other original material on the early 
work of the Wright brothers. His donation made our 
Archives the most complete and unique depositary of 
Wright aeronautica in existence. 

‘‘When the Institute first considered the inauguration 
of an annual lecture here, comparable to the Wilbur 
Wright Memorial Lecture presented each year by the 
Royal Aeronautical Society in England, we naturally 
conferred with Mr. Orville Wright to determine his 
wishes and opinion on the subject. He expressed his 
approval of our plan to establish the Wright Brothers 
Lecture and attended the first of these lectures which 
was given here on December 17, 1937, by Professor B. 
Melvill Jones, another of our distinguished English 
members. As you all know, Mr. Wright very seldom 
makes public appearances or statements of any kind, 
particularly statements as to his opinions on any public 
matters—aeronautical or general. Therefore, I feel 


that the Institute is signally honored to have received 
from Mr. Wright the following telegram sent to us on 
this occasion commemorating the thirty-eighth anni- 
versary of the first flights at Kitty Hawk: 


“Please extend to Professor Southwell my cordial 
greetings and an expression of my belief that the aero- 
plane which, in the hands of barbarous men, has been 
doing England and the world inestimable harm will yet 
be instrumental in establishing for the entire world a 
guarantee of human liberty and of lasting peace. 

ORVILLE WRIGHT 


“T now have another announcement which I feel is 
appropriate to make before this gathering of members 
and friends of the Institute who are all greatly interested 
in every forward step made by our organization. Two 
years ago, Mr. Paul Kollsman, a Fellow of the Institute, 
and the Square D Company of Detroit, gave us $50,000 
to be spent on the establishment and maintenance of an 
aeronautical lending library. The Paul Kollsman 
Library has rendered such effective aid in lending books 
by mail to aeronautical students and technical special- 
ists all over the country that Mr. Kollsman has decided 
to establish this service on a permanent basis. He in- 
formed me this week that he will make an additional 
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gift to the Institute so as to establish a permanent fund 
of $65,000 to be used for the maintenance of the library. 
I am sure that all of you here join with me in expressing 
to Mr. Kollsman the appreciation of the Institute for his 
great generosity in endowing an activity which is so 
valuable to American aeronautics in these critical 
times.”’ 

Dr. Hunsaker: “We have the honor to have 
present the official representative of British science 
and of her learned societies, Dr. Charles Darwin, 
Chief of the British Central Scientific Office in Wash- 
ington. Dr. Darwin implements the close collabora- 
tion and exchange of ideas between research workers 
of this country and those of the British Empire, which 
is proving so effective in the creation of new weapons 
and techniques.”’ 

Dr. Darwin: ‘‘I hardly need to tell you my pleasure 
at being able to be with you tonight and the distinction 
that I feel has been conferred upon me by the ap- 
pointment as Honorary Chairman of this Wright 
Brothers Lecture. When I came to this country nine 
months ago in connection with my work in Washing- 
ton, I came as a friend. When my duties take me 
back to England, as they will in the near future, I am 
glad that I shall leave you not only as a friend, but as 
an ally.” 

Dr. Hunsaker: ‘‘Now a word of introduction for 
our guest, Professor Southwell. You already know 
he is important to us because he was invited to cross 
the Atlantic to speak to us. You can also realize this 
winter journey was not undertaken in response to a 
mere social invitation. The authorities, including 
Dr. Darwin, who control such journeys made the visit 
possible. 

“Professor Southwell brings to any problem a broad 
experience in practically all of the aeronautical sci- 
ences. He illustrates the ideal definition of a Fellow 
of this Institute. He has been for many years a mem- 
ber of the British Aeronautical Research Committee 
and more recently has been directing the work of the 
Civil Defense Research Committee. His many con- 
tributions to the development of engineering science 
have been characterized by the approach of the ana- 
lyst with fluent command of the known mathematical 
tools, and at the same time with a clear recognition of 
physical realities 

“The Institute welcomes Professor Southwell as its 
1941 Wright Brothers Lecturer, and I am honored to 
present him to you now.” 

Professor Southwell expressed his thanks to the 
Chairman and to the Institute for the opportunity to 
come to the United States to give his paper and re- 
ferred to previous pleasant visits here. He then pre- 
sented the Fifth Wright Brothers Lecture. 
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OFFICERS ELECTED FOR 1942 
Hatt L. HIBBARD TO BE PRESIDENT 


The Council has elected Hall L. Hibbard, Vice President and 
Chief Engineer of Lockheed Aircraft Corporation, as President of 
the Institute for 1942. He will be inducted into office at the 
Honors Night Dinner held in New York on January 27. 

Mr. Hibbard is responsible for the design of a long series of 
the noted Lockheed airplanes. He helped in pioneering the 
development of the twin tail and the single-spar all-metal airplane 
wing. His designs have used increasingly heavy wing loadings 
and include an experimental tail-first or canard type of airplane. 
In the development of air transport types and the Army’s famous 
stratosphere airplane, his work has been outstanding. 

Hall L. Hibbard was born in Fredonia, Kansas, and spent his 
boyhood in the Philippine Islands, where his father was a mis- 
sionary for 41 years. After graduation from the College of Em- 
poria, he took a postgraduate course in aeronautics at the Massa- 
chusetts Institute of Technology and received an M.A. degree in 
aeronautical engineering in 1927. He became Chief Engineer of 
the Lockheed company in 1932, having previously been on the 
engineering staff of Stearman Aircraft Company since his gradua- 


tion. 
Other officers elected by the Council are as follows: 


VICcE-PRESIDENTS 


C. S. Jones, President, Casey Jones School of Aeronautics 
Grover Loening, Consulting Aeronautical Engineer 
J. G. Vincent, Vice-President of Engineering, Packard Motor 


Car Company 
Thomas A. Morgan, President and Director, The Sperry Cor- 


poration 
EXECUTIVE VICE-PRESIDENT—Lester D. Gardner 
SECRETARY—C. E. Sinclair 
TREASURER—Earl D. Osborn, President, Edo Aircraft Corpora- 
tion 
ASSISTANT TREASURER—George R. Forman 
CONTROLLER— Joseph J. Maitan 


FELLOwS FoR 1941 ELECTED 


The following members were elected Fellows of the Institute in 
1941: 
George W. Brady, Chief Engineer, Curtiss Propeller Div., 
Curtiss-Wright Corp. 

John W. Crowley, Jr., Sr. Aero. Engineer, Langley Mem. Aero. 
Lab., N.A.C.A. 

Ralph S. Damon, Pres., Republic Aviation Corp. 

William K. Ebel, Vice-Pres. in charge of Engineering and Chief 
Engineer, The Glenn L. Martin Co. 

Michael E. Gluhareff, Chief of Design and Aerodynamics, 
Vought-Sikorsky Aircraft Div., United Aircraft Corp. 

R. M. Hazen, Vice-Pres. and Chief Engineer, Allison Div., 
General Motors Corp. 

Comdr. ¥. W. Pennoyer, Jr., U.S.N., Head of Engineering 
Div., Bureau of Aeronautics, Navy Department. 

Harry A. Sutton, Aero. Engineer, Consolidated Aircraft Corp. 

Capt. W. W. Webster, U.S.N., Manager, Naval Aircraft Fac- 
tory, Philadelphia. 

Robert J. Woods, Chief Engineer, Bell Aircraft Corp. 

Qualifications required of a Fellow of the Institute are that he 
shall have attained a position of distinction in aeronautics and 
have made notable and valuable contributions in one of the 


aeronautical sciences or aeronautical engineering. Not more 


than ten members of the Institute may be elevated to this grade 
of membership each year. 

Selection of those to receive this honor is made by the present 
American Fellows and Honorary Fellows of the Institute. 
Nominations are submitted by them and Fellows selected from 
this list by secret elimination ballots. 


NEW CORPORATE MEMBERS 


The Institute is pleased to announce the affiliation of two 
nationally known financial institutions which became Corporate 
Members at the beginning of 1942: 

The Chase National Bank of the City of New York, which has 
for some time been cooperating with all branches of the aviation 
industry and airline operators. It has made loans for the ex- 
pansion of plant facilities, working capital needs, and for the ac- 
quisition of transport equipment. 

The National City Bank of New York, which has long been 
actively identified with the aircraft industry. 

The Institute appreciates the support and interest of these 
organizations in its work and will be able to assist them, as it does 
all its Corporate Members, to keep in closer touch with the entire 
aviation industry and its rapidly changing technical aspects. 


KOLLSMAN FUND MADE PERMANENT 


At the Wright Brothers Lecture on December 17 it was an- 
nounced that Paul Kollsman, Fellow and Benefactor Member of 
the Institute, had made an additional gift to the Institute to set 
up a permanent endowment fund for the support of The Paul 


Kollsman Library. 





Paul Kollsman 


In 1940 Mr. Kollsman and the Square D Company of Detroit 
gave $50,000 to the Institute to establish a library to lend books 
to members and others wishing to increase their knowledge in 
any of the aeronautical sciences. Now Mr. Kollsman has given 
the Institute an additional $25,000, which, with the remainder 
of the original fund, will provide a permanent endowment fund of 
$65,000 to support the activities of The Paul Kollsman Library. 

The widespread use of the library since its establishment, 
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through the hundreds of aeronautical books that have been 
loaned by mail without charge to aeronautical students and engi- 
neers throughout the country, attests to the value of the con- 
tribution that Mr. Kollsman’s generous gifts are making to the 
country’s air defense. His great contributions through the 
widely used aircraft instruments which he has designed and de- 
veloped during the past fifteen years are well known to every one 
in aviation. The opportunities for the dissemination of aero- 
nautical knowledge which his donation to the Institute has made 
possible are likewise appreciated by the aeronautical profession 
and will be long remembered. 


DONATION FROM R. H. FLEET 


Major R. H. Fleet, of San Diego, California, founder of the 
Consolidated Aircraft Corporation, and a Founder member of the 
Institute, made a gift of $10,000 to the Institute last month. 
+: Inthe telegram in which Major Fleet notified the Institute of his 
intention to make this donation, he said: ‘‘Being deeply inter- 
ested in aeronautical science and realizing the fine work your Insti- 
tule is carrying on along these lines, J desire to give $10,000 to assist 


og 





R. H. Fleet 


The Officers and Council at their last meeting passed a resolu- 
tion expressing, on behalf of the Institute, their appreciation of 
Major Fleet’s interest in the progress of the Institute as stated in 
his message and so generously demonstrated through his gift for 
the furtherance of its work. 


FAIRCHILD PUBLICATION FUND ESTABLISHED 


Sherman M. Fairchild, Founder member of the Institute, 
Chairman of the Board of the Fairchild Aviation Corporation 
and the Fairchild Engine and Airplane Corporation, donated to 
the Institute last month a fund of $25,000 for the publication 
of technical papers, reports and other aeronautical material that 
would not otherwise be published in aviation trade magazines 
or scientific journals. 

In its resolution accepting this gift and conveying the apprecia- 
tion of the Institute to Mr. Fairchild for his generosity, the Coun- 
cil set it up as a fund to be known as ‘‘The Sherman Fairchild 
Publication Fund’’ and stipulated that all publications issued by 
it should bear this imprint. The fund will be managed by a com- 
mittee consisting of Mr. Fairchild, Charles H. Colvin and Lester 
D. Gardner. They will be advised as to the choice of material 


for publication by the Editorial Board of the Journal of the Aero- 
nautical Sciences. 

It is intended that part of the fund be used in the near future 
to publish certain technical reports whose value, though vital to 
the aeronautical sciences, is not of sufficient general interest to 
warrant its inclusion in magazines of general circulation. It will 
also make possible the reproduction of engineering papers whose 
circulation, for reasons of military secrecy, must be restricted to 
the government approved list of technical people engaged in 
national defense work who are authorized to receive confidential 
material. The services of the fund have already been offered to 
the Army and Navy air forces for the publication and distribution 
of reports as they may suggest. 


Mayor Harry G. ARMSTRONG 
NAMED FOR JEFFRIES AWARD 


Major Harry G. Armstrong, U.S. Army Medical Corps, in 
charge of research at the School of Aviation Medicine, Randolph 
Field, Texas, has been named to receive the John Jeffries Award 
given by the Institute, for outstanding contributions to the 
advancement of aeronautics through medical research. 

Major Armstrong established the Aero Medical Research 
Laboratory of the Army Air Corps at Wright Field, Dayton, 
Ohio, in 1934 and was the director of this unit until his recent 
assignment to duty at Randolph Field. Born in DeSmet, 
South Dakota, February 17, 1899, he took pre-medical courses 
at the University of Minnesota, received a B.S. degree from the 
University of South Dakota in 1923 and his M.D. from the 
University of Louisville (Kentucky) in 1925. He served as an 
interne at St. Luke’s Hospital, St. Paul, Minnesota, and is a 
graduate of the Army Medical School, the Army Medical Field 
Service School and the Army School of Aviation Medicine. 

The John Jeffries Award will be conferred on Major Armstrong 
in recognition of his pioneering studies done on the physiological 
and psychological effects of flying at high altitude and in high 
speed maneuvers. He was one of the first to recognize and ac- 
curately describe some specific medical results of flying, such as 
aero-neurosis, a kind of mental and physical fatigue experienced 
by flyers under certain conditions; aero-otitis, an effect of high 
altitude and acrobatic flying on the human middle ear; and aero- 
embolism, a reaction similar to the ‘“‘bends’’ suffered by deep sea 
divers which affects pilots who climb to high altitudes too 
rapidly and without proper safeguards. 

Most of his researches on the physiological effects of low baro- 





Sherman M. Fairchild 
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Honors Night Dinner 
Horet Watporr-Astoria, New York 
TueEspDAY, JANUARY 27, 7:30 P.M. 


Marking the Institute’s tenth year as the only scientific aero- 
nautical society in the United States, the Honors Night Dinner 
opens the sessions of the Tenth Annual Meeting in New York. 
Due to the large attendance at previous dinners and the even 
greater number expected on this anniversary occasion, more am- 
ple accommodations have been engaged for this affair than in 
previous years. Reservations for the dinner are $6.00 per person. 
The program for the Dinner follows: 
Greeting to Members and Guests, by Frank W. Caldwell, 
President of the Institute of the Aeronautical Sciences. 
Presentation of The Daniel Guggenheim Medal to Juan T. 
Trippe. 

Presentation of The Sylvanus Albert Reed Award to Dr 
Theodore von Karman. 

Presentation of The Octave Chanute Award to Melvin N. 
Gough. 

Presentation of The Lawrence Sperry Award to Ernest G. 
Stout. 

Presentation of The John Jeffries Award to Major Harry G. 
Armstrong (M.C.). 

Presentation of The Robert M. Losey Award (recipient to be 

announced). 

Presentation of Certificates to Fellows of the Institute Elected 

in 1941. 

Presentation of Certificates of Honorary Membership in the 

Institute to 

Hon. Robert A. Lovett, Assistant Secretary of War for Air; 

Hon. Artemus L. Gates, Assistant Secretary of the Navy for 
Air: 

Hon. Robert H. Hinckley, Assistant Secretary of Commerce 
for Air. 

Induction of Hall L. Hibbard as President of the Institute for 
1942. 


Technical Sessions of Annual Meeting 
Pupin Puysics LABORATORIES 
CotumsiA UNIversity, New York 


Technical papers will be presented at morning, afternoon and 
evening sessions over a period of four days. Luncheon will be 


metric pressure, low temperature and lack of oxygen were carried 
out in the high altitude test chamber of the Aero Medical Re- 
search Unit at Wright Field. Major Armstrong customarily 
made himself the first human subject for a new test, often under 
experimental conditions which were dangerous to his own life 
His work has done much to stimulate further research in the 
medical aspects of aviation and has led to the development of 
oxygen supply and pressurizing apparatus and other precaution- 
ary devices and procedures for safeguarding pilots against the 
physical effects of modern military flying. Numerous reports on 
his work have been published in medical and aeronautical 
journals. His book on “Principles and Practice of Aviation 
Medicine” is the most complete and authoritative text yet pub- 
lished on the subject. A 

Major Armstrong is a member of the Aero Medical Association, 
the Association of Military Surgeons and the American Medical 
Society. He is an Associate Fellow of the American College of 
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TENTH ANNUAL MEETING AND Honors NIGHT DINNER 


JANUARY 27-31, 1942, NEW YORK 








SCIENCES 







served each day at the Columbia University Men’s Faculty Club. 
Since numerous changes have had to be made in the list of papers 
to be presented, the titles are not given here. Members of the 
Institute have been informed as to the papers being presented 
through programs sent by mail. Final program will be available 
at the meeting. The schedule and location of sessions is given 
below. 

WEDNESDAY, JANUARY 28 
Structures, Part 1, Room 401 
Chairman, Eugene E. Lundquist 
Power Plants, Room 428 
Chairman, R. N. DuBois 
Structures, Part 2, Room 401 
Chairman, John E. Younger 
Radio and Instruments, Room 428 
Chairman, W. A. Reichel 
Annual Business Meeting, Room 401 


9:15 a.m. 


bo 


:00 p.m. 


or 


:00 p.m. 


oo 


Design, Room 401 
Chairman, Fred E. Weick 


:00 p.m. 


Tuurspay, JANUARY 29 


Materials, Part 1, Room 401 

Chairman, C. F. Nagel, Jr. 

Rotating Wing Aircraft, Room 428 
Chairman, Ralph H. McClarren 
Materials, Part 2, Room 401 

Chairman, A. V. de Forest 

Physiologic Problems, Part 1, Room 428 
Chairman, Arnold D. Tuttle 

Physiologic Problems, Part 2, Room 401 
Chairman, John G. Jenkins 


9:15 a.m. 


bo 


:00 p.m. 


o 


:00 p.m. 


Fripay, JANUARY 30 


Aerodynamics, Part 1, Room 428 
Chairman, A. M. Kuethe 
Meteorology, Part 1, Room 401 
Chairman, C. F. Brooks 
Aerodynamics, Part 2, Room 428 
Chairman, William Bollay 
Meteorology, Part 2, Room 401 
Chairman, F. W. Reichelderfer 


:15 a.m. 


:00 p.m. 


bo 


8:00 p.m. Air Transport, Room 401 


Chairman, Charles Froesch 


Physicians and of the Institute of the Aeronautical Sciences. 
He was chosen to receive the John Jeffries Award for 1941 by 
a committee consisting of Dr. G. W. Lewis, Director of Research 
of the National Advisory Committee for Aeronautics; Gill Robb 
Wilson, President of the National Aeronautical Association; 
Edgar S. Gorrell, President of the Air Transport Association 
of America; Capt. J. R. Poppen, U.S.N., President of the Aero 
Medical Association; Frank W. Caldwell, President of the 
Institute of the Aeronautical Sciences; and Lester D. Gardner, 
President of the Aeronautical Archives of the Institute. 

The award is named after Dr. John Jeffries, a Boston physician, 
who was the first American to make scientific observations from 
the air. He made a flight from London with Blanchard, the 
French balloonist, in 1784 and took recordings of the tempera- 
ture, humidity and density of the air. In January, 1785, Dr. 
Jeffries and Blanchard made the first crossing of the English 
Channel by air. 
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LAWRENCE SPERRY AWARD FOR 1941 TO 
ERNEST G. STOUT 


The Lawrence Sperry Award for a notable contribution made 
by a young man to the aeronautical sciences will be presented to 
Ernest G. Stout at the Institute’s Honors Night Dinner in New 
York on January 27. The award consists of an honorarium of 
$250 and a certificate which will bear the citation ‘‘For his con- 
tributions to the experimental determination of the hydro- 
dynamic stability of model flying boats and seaplanes.”’ 

Mr. Stout, a MEMBER of the Institute, is Engineer in charge 
of Aerodynamics and Flight Test for Consolidated Aircraft Cor- 
poration. He was an Aerodynamic Engineer for this company 
from 1936 until 1938, then Engineer in charge of Hydrodynamics 
until he recently assumed his present position. 

Born in Moran, Kansas, April 25, 1913, he graduated from 
Western High School, Washington, D. C., in 1931, and received 
The Perry Collins McDonough Scholarship for four years at the 
Daniel Guggenheim School of Aeronautics, New York University. 
He received the degree of B.S. in Mechanical Engineering from 
New York University in 1935 and served as a Research Assistant 
there for three months in the summer of that year. In 1939 Mr. 
Stout received a non-resident degree as Aeronautical Engineer 
for which he submitted a thesis entitled “‘Charts for the Deter- 
mination of Seaplane Performance.” 

He was an Aerodynamic Engineer with The Glenn L. Martin 
Company from July, 1935, until November, 1936, when he joined 
the engineering staff of Consolidated. 

Mr. Stout’s work on the hydrodynamic stabtlity of flying boats 
and seaplanes, for which he receives the award, was begun in 1938. 
He read a general paper on ‘‘Notes on the Design and Stability of 
Flying Boats’’ before the Los Angeles Section of the Institute in 
September, 1939, and his report on ‘‘Experimental Determination 
of Hydrodynamic Stability’’ was published in the Journal of 
December, 1940. These described his design of seaplane and fly- 
ing boat models with controllable elevator angle and certain 
modifications on wing and tail surfaces to counteract scale effect. 
The methods of mounting and testing these models as supported 
from a boat in open water or from the carriage in a towing basin, 
have permitted the determination of the effect on stability of the 
center of gravity location and hull bottom modifications. Thus 
the changes in hull design and construction necessary to prevent 
dangerous porpoising can be made and tested in small models be- 
fore being incorporated in a full size flying boat. Mr. Stout’s 
methods have been used by three government departments and 
several manufacturers and have been closely checked by full 
scale tests. 

The Award, honoring the memory of Lawrence Sperry, pioneer 
aviator and inventor who died in 1923 at the age of 31, is endowed 
with a fund of $10,000 by his brothers and sister. 

The recipient of The Lawrence Sperry Award is chosen by a 
Board of Award whose members in 1941 were Lt. Col. James H. 
Doolittle, Lester D. Gardner, Charles L. Lawrance, Grover 
Loening, Glenn L. Martin and Elmer A. Sperry, Jr. 


News oF INSTITUTE MEMBERS 


Eli Amaneul, Technical Member, has joined the staff of the 
Munk Aeronautical Laboratory as a member of the firm. 

Albert A. Andrake, A.F.I.Ae.S., Project Tool Engineer, The 
Glenn L. Martin Co., for the past two years, has been transferred 
to The Glenn L. Martin-Nebraska Co., Omaha, Nebraska, as 
Chief Production Engineer. 

W. John Beil, Technical Member, is now in the Aerodynamics 
Section, Airplane Division, Curtiss-Wright Corp., Buffalo, N. Y. 

Raymond L. Bisplinghoff, Technical Member, formerly Re- 
search Assistant, has become Instructor in Aeronautical Engi- 
neering at the University of Cincinnati. 

Allan F. Bonnalie, A.F.I.Ae.S., Assistant to United Air Lines’ 
Executive Vice-President in charge of operations, has taken a 


one year leave of absence effective since Dec. 1, to join the 
Bureau of Aeronautics, Navy Dept. He is a Lt. Comdr. in the 
U.S. Naval Reserve. 

Douglas Campbell, M.I.Ae.S., has been named Vice-President 
of Pan American-Grace Airways. He has been active in Panagra 
since May, 1935. 

Robert Alden Cornog, Technical Member, who has been a 
Research Fellow at the Univ. of California, is now a Physicist, 
Bureau of Ordnance, Navy Dept. 

Peter de Florez, Technical Member, has left the Ranger Air- 
craft Engines Div. of Fairchild Engines & Airplane Corp. to 
join the Link Aviation Devices at Binghamton, N. Y. 

Albert B. Deyarmond, M.I.Ae.S., is now a Captain in the 
Army, Specialist Reserve U.S. Army Air Corps, at Wright Field. 
He was formerly Chief of the Stress Dept., Vultee Aircraft, Inc. 

Lt. Col. James H. Doolittle, F.I.Ae.S. and Past President of 
the Institute, has been promoted to this rank recently, it was an- 
nounced by the War Department. Formerly Director of Avia- 
tion for Shell Oil Company, he has been on active duty since July, 
1940, as a Major in the Procurement Division of the Army Air 
Corps. From November of that year he was stationed at Detroit 
and is now engaged on research and engineering problems per- 
taining to aircraft production, with headquarters at Baltimore. 

Ralph H. Draut, M.I.Ae.S., has been named Chief of Struc- 
tures, The Glenn L. Martin-Nebraska Co. Formerly a Project 
Stress Engineer for The Glenn L. Martin Co., he has been as- 
sociated with such firms as the Curtiss Aeroplane & Motor Co., 
Stearman Aircraft Co., and Waco Aircraft Co. 

John C. Dykes, Technical Member, formerly the Assistant to 
the Research Engineer of British Overseas Airways Corp. at 
Bristol, England, has been transferred to Asmara, Italian East 
Africa (British Occupied). He is a Technical Officer at the 
company’s base there. 

Reuben H. Fleet, M.I.Ae.S., having sold his interest in the 
Consolidated Aircraft Corp., which he founded, to Vultee Air- 
craft, Inc., has resigned as President and Director of the com- 
pany. He will remain with Consolidated in an advisory ca- 
pacity. 

A. Milton Drummer, Technical Member, is now a Stress 
Analyst with the Beech Aircraft Corp. 

George L. Fruauff, Jr., Technical Member, recently with the 
Brewster Aeronautical Corp., has become associated with the 
Western Electric Co. at Kearney, N. J. 

R. E. Gillmor, A.F.I.Ae.S., President of the Sperry Gyro- 
scope Co., has become a member of the Board of Directors of the 
Fairchild Engine and Airplane Corp. 

George M. Goldman, Technical Member, has joined the 
Materiel Division, U.S. Army Air Corps, Wright Field. 

Nicholas J. Hoff, A.F.I.Ae.S., Assoc. Professor of Aeronautical 
Engineering at the Polytechnic Institute of Brooklyn, has been 
elected an Associate Fellow of the Royal Aeronautical Society. 

Roger Wolfe Kahn, M.I.Ae.S. and member of the Council, 
has been appointed a member of the Aviation Planning Staff, a 
research and advisory body attached to the Civil Air Patrol. 

Major Reed G. Landis, M.I.Ae.S., of the Office of Civilian De- 
fense, which organized the Civil Air Patrol, is head of the Plan- 
ning Staff. 

Herbert F. Kueck, Technical Member, formerly with the 
Beech Aircraft Corp., is now an Instructor at the Guggenheim 
School of Aeronautics, N. Y. U. 

Harvey A. Leich, Technical Member, recently with the Grum- 
man Aircraft Engineering Corp., is now a 2nd Lt. and Engi- 
neering Officer with the First Marine Aircraft Wing, U.S.M.C. 
Reserve. 

Richard C. du Pont, M.I.Ae.S., President of All American 
Aviation, Inc., of Wilmington, Del., has announced that Grover 
Loening, F.I.Ae.S., aircraft manufacturer and engineer, has 
become associated with the company as Engineering Advisor, and 
also has been elected a member of the Board of Directors. 
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Cameron N. Lusty, Technical Member, has been named Chief 
Engineer of the Mercury Aircraft Co. He was formerly a 
Stress Analyst at the Taylorcraft Aviation Corp. 

Victor J. Martin, M.I.Ae.S., has joined Vultee Aircraft, Inc., 
as an Aerodynamics Research Engineer. For the past four 
years he has been teaching aeronautical engineering at the 
University of Washington and at the time of his resignation was 
an Asst. Professor. 

Robert A. McConnell, Technical Member, is now a Research 
Associate, Radiation Laboratory, Mass. Inst. of Technology. 

Harold P. Moon, M.I.Ae.S., formerly General Manager, 
Summit Aeronautical Corp., is now a test pilot of military air- 
craft, Airplane Div., Curtiss-Wright Corp., Buffalo, N. Y. 

Robert R. Osborn, F.I.Ae.S., is now Asst. Engineering Ad- 
visor of the Bell Aircraft Corp. A graduate of the Univ. of 
Pennsylvania, Mr. Osborn has been with the Curtiss Aeroplane 
& Motor Co., and more recently was Asst. Director of Engineer- 
ing, Vultee Aircraft, Inc. 

Kingsbury E. Parker, Jr., Technical Member, has been made 
a Captain in the Royal Air Force Ferry Command, Montreal. 

C. Kendall Perkins, A.F.I.Ae.S., has left American Airlines 
and is now a Project Engineer with the McDonnell Aircraft 

«Corp. 

Capt. John R. Poppen, U.S.N., M.I.Ae.S., has been transferred 
from the Bureau of Aeronautics, to the Bureau of Medicine 
and Surgery, Navy Dept. 

Allen E. Puckett, Technical Member, is now a Research As- 
sistant, Calif. Inst. of Technology. 

E. Robert Reff, M.I.Ae.S., formerly Structural Engineer at 
the Boeing Aircraft Co., is now Structural Engineer with the 
Aircraft Div. of Fairchild Aircraft Corp. at Hagerstown, Md. 

Peter F. Rossmann, M.I.Ae.S., formerly Development and 
Research Engineer, has been named Assistant to the Director of 
Military Engineering, Airplane Div., Curtiss-Wright Corp., 
Buffalo, N. Y. 

Alexander Satin, M.I.Ae.S., has joined the Consolidated Air- 
craft Corp. as Design Engineer. He was formerly Chief Engi- 
neering Instructor, School of Aeronautics, California Flyers, 
Inc. 

Richard A. Schram, Technical Member, has been named 
Senior Ground School Supervisor, CAA, in the Chicago area. 

Comdr. C. H. Schildhauer, U.S.N., A.F.I.Ae.S., recently 
Marine Equipment Sales Manager for The Glenn L. Martin Co., 
reported for active duty with the Navy on December 1. 

N. F. Silsbee, M.I.Ae.S., aviation writer and lecturer, and 
formerly Comptroller of Beaver College, Jenkintown, Pa., has 
been on active duty since October as a Captain in the Air Corps. 
He is Chief of the Research Section, Public Relations Branch, 
Hq. Army Air Forces. 

Lt. Emil L. Sorenson, Technical Member, is now a Test Pilot, 
Materiel Div., U.S. Army Air Corps, Wright Field. 

Morril B. Spaulding, Technical Member, has joined The Glenn 
L. Martin Co. as a Project Liaison Engineer. He was formerly 
with the Pan American Airways System as a Jr. Engineer. 

John F. Vogel, Technical Member, has been called to active 
duty as a Ist Lieut., Signal Corps, MacDill Field, Florida. He 
was formerly with the Stinson Aircraft Div. of Vultee Aircraft, 
Inc. 

Capt. W. W. Webster, U.S.N., F.I.Ae.S., has recently been 
transferred from the office of the Material Branch of the Navy 
Bureau of Aeronautics in Washington to become Manager of the 
Naval Aircraft Factory in Philadelphia. 

Robert B. White, M.I.Ae.S., formerly Asst. to Vice-Pres., 
TWA, Inc., is now Factory Manager, North American Aviation, 
Inc., of Kansas. 

Robert S. Wolfsohn, Technical Member, now assigned to a 
project at the National Bureau of Standards in Washington, is 

an Assistant Mechanical Engineer on the staff of the N.A.C.A. 
Bradley H. Young, Technical Member, who has been connected 
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with the Maintenance staff of Pan American-Grace Airways 
since 1937, has been appointed Maintenance Engineer at Lima, 
Peru. 

Shao W. Yuan, Technical Member, has joined the Vidal 
Research Corp. as an Aeronautical Engineer. 

Harry Zuckerberg, M.I.Ae.S., is now Chief of Research and 
Development, Platt-LePage Aircraft Co. He was formerly 
Chief Structures Engineer. 


PERSONNEL OPPORTUNITIES 


The Personnel Bureau serves individual members, as well as 
organizations seeking to employ aeronautical specialists. Any 
member or organization may have requirements listed without 
charge. 


Available 


Aeronautical and Mining Engineer, familiar with Air Corps 
and Bureau of Aeronautics contract procedure desires an op- 
portunity to make new connection. Continuously employed for 
approximately 15 years on production engineering work involv- 
ing production schedules, materials and other aspects of pro- 
curement planning incident to aircraft production. Address 
reply to Box 146, Institute of the Aeronautical Sciences. 


STUDENT BRANCHES 
Aeronautical University. A paper on ‘‘The Workings and 
Uses of the Terrain Clearance Indicator and the Flightray’”’ was 
presented by Edward Bedrosion at a meeting held on December 2. 

Indiana Technical College. Two films, borrowed from the 
Institute, were shown at a meeting held on December 2. They 
were “Conquest of the Air’? and ‘‘Goodrich De-Icer Develop- 
ment.”’ 

Massachusetts Institute of Technology. W. A. Petrasek, of 
American Airlines, gave a talk on the Sperry “Flightray”’ at a 
meeting held on December 9. At a previous meeting, Thomas 
D. Perry, of The Resinous Products & Chemical Company, gave 
a talk and demonstration of aircraft plywoods. 

University of Pittsburgh. At the first meeting of this new 
Student Branch, officers were elected as follows: Marshal 
Nesbitt, Joseph Siro, Vice-Chairman; Raymond 
Zavasky, Treasurer; Theodore Just, Secretary. J. Elwood 
Arnold is Faculty Advisor. 

Rensselaer Polytechnic Institute. At a meeting held on De- 
cember 11, two U.S. Army Air Corps films were shown: ‘‘Train- 
ing of Air Corps Cadets” and ‘Wings of Steel.’”’ Later in the 
evening, the present officers of the Branch were re-elected. 
They are: George McTigue, Chairman; Warren Tucker, Treas- 
urer; and Richard P. Graham, Secretary. 


Chairman; 


CHANGES IN MEMBERSHIP 


The following changes in the membership of the Institute 
have occurred since the publication of the previous issue of the 


Journal. 
TRANSFERRED TO ASSOCIATE FELLOW GRADE 


Blount, Earl Ellsworth, Chief of Structures, Airplane Div., 
(St. Louis Plant), Curtiss-Wright Corp. 

Doman, Carl Templar, B.S. in E.E.; Vice-Pres. in charge of 
Engineering, Aircooled Motors Corp. 

Hoff, Nicholas John, Dipl. Ing.; Assistant Professor of Aero 
Eng., Polytechnic Institute of Brooklyn. 


ELECTED TO MEMBER GRADE 


Barton, Millard Vernon, Ph.D.; Asst. Prof. of Mech. Engi- 
neering, University of Maryland. 
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Goodwin, Ralph Talbot, Ph.D.; Manager, Aviation Dept., 
Shell Oil Co., Inc. 

Greer, Edward Maurice, B.S. in Ae.E.; Chief Design Engineer, 
Simmonds Aerocessories, Inc. 

Horrisberger, John Frederick, B.S. in M.E.; Project Engi- 
neer (Airplanes), U.S. Naval Aircraft Factory, Philadelphia. 

Howard, Henry Bernon, B.A. in M.E.; B.Sc.; Assistant Direc- 
tor, British Air Commission. 

Kelly, Everett Lowell, Ph.D.; Associate Professor of Psy- 
chology, Purdue University (See A. M. of S.). 

Krzywoblocki, Maria Zbigniew, Dipl. Ing.; M.S.; A.F.R.Ae.- 
S., Asst. to Supervisor of Tool Control & Works Method Dept., 
Aircraft Div., Canadian Car and Foundry Co. 

Shaw, William Ulric, B.A.Sc.; Chief Engineer, Aircraft Div., 
National Steel Car Corp., Ltd., Canada. 

Watson, Donald Reed, Chief Weight Engineer, Fleetwings, 
Inc. 

Wiggins, Garwood Owen, Executive Engineer, Ranger Air- 
craft Engines Div., Fairchild Engine & Airplane Corp. 


TRANSFERRED TO MEMBER GRADE 


Allen, Fred Charles, M.E. in Ae.; Ass’t. Chief of Stress, El 
Segundo Div., Douglas Aircraft Co. 

Walker, Ernest E., A.B.; M.L.A.; Aeronautic Constructor for 
Airports and Airways. 


ELECTED TO INDUSTRIAL GRADE 


Zweng, Charles Alfonso, LL.B.; Owner and Manager, Pan 
American Navigation Service. 


ELECTED TO TECHNICAL MEMBER GRADE 


Anderson, Benjamin Monroe, B.A.; Detail Engineer, Boeing 
Aircraft Co. 

Bates, William Elmer, Jr., B.S. in M.E.; Technical Engineer- 
ing Assistant, Blade Design, Curtiss Propeller Div., Curtiss- 
Wright Corp. 

Ebeling, Wilbur Ludlow, B.S. in Ae.E.; Instructor, Wayne 
University. 

Greenewald, Raymond Eugene, B.S. in Ae.E.; Instructor 
Aero. Dept., Tri-State College. 

Hermes, Raymond Francis, B.Ae.E.; Stress Analyst, Aeronca 
Aircraft Corp. 

Ill, Edward Joseph, II, Weight and Balance Engineer, AGA 
Aviation Corp. 

Lentz, George Lawrence, B.S. in E.E.; Engineer, General 
Electric Co. 

Meyer, Carl George, Engineer, Republic Aircraft Corp. 

Neeley, Glen E., Assistant Instructor, Post Schools Dept., 
U.S. Army Air Corps, McClellan Field, Calif. 

Pfeiffer, Shelly Bernhardt, B.S. in M.E.; Stress Group, El 
Segundo Div., Douglas Aircraft Co. 

Schnaubelt, Franz Joseph, Planning Engineer, Group Leader, 
S. I. Group, Production Engineering Dept., Consolidated Air- 
craft Corp. 

Schell, Charles Joseph, Jr., Engineering Instructor, Casey 
Jones School of Aeronautics. 

Westergaard, John William, Engineer, Boeing Aircraft Co. 


TRANSFERRED FROM STUDENT MEMBER TO TECHNICAL MEMBER 
GRADE 


Alexander, Francis Brown, B.S. in M.E.; 2nd Lt. A.C., Air 
Corps Tactical School. 

Belzer, Maurice, B.S. in Ae.E.; Engineer, Vought-Sikorsky 
Aircraft Div., United Aircraft Corp. 

Bemer, Robert W., A.B.; Engineering Calculator, Douglas 
Aircraft Co. 


Beven, Kenneth A., Engineering Draftsman, El Segundo Div., 
Douglas Aircraft Co. 
Bielat, Ralph Philip, B.S. in Ae.E.; Jr. Aero. Engineer, Lang- 
ley Memorial Lab., N.A.C.A. 
Blake, Robert W., S.B. in Ae.E.; Member Engineering Dept., 
Pan American Airways, Inc. 
Bogert, Robert C., B.S. in Ae.E.; Aviation Cadet, U.S. Army 
Air Corps. 
Bonney, Robert Wm., Inspector, Consolidated Aircraft Corp. 
Brogger, Charles, C.B.S. in Ae.E.; Draftsman, Designer, Ex- 
perimenter, Camfield Mfg. Co. 
Buie, John C., Instrument Technician, Mackenzie Air Service 
Ltd. 
Burgeson, John Edward, Draftsman, Ryan Aeronautical Co. 
Castiglione, Frank P., B.S. in Ae.E.; Engineering Draftsman, 
Dept. of Commerce. 
Chin, Charles, B.S. in Ae.E.; Loftsman, Airplane Div. (Co- 
lumbus Plant), Curtiss-Wright Corp. 
Christie, Sidney A., Draftsman, Ryan Aeronautical Co. 
Daughaday, Hamilton, M.S. in Ae.E.; Structures Dept., 
Airplane Div. (Buffalo Plant), Curtiss-Wright Corp. 
Dilworth, John A., B.S. in Ae.E.; Lt. U.S. Army Air Corps. 
Dugan, James C., A.E.; Jr. Aero. Engineer, U.S. Army Air 
Corps. 
Frei, Ralph E., Aero. Engineer, The Glenn L. Martin Co. 
Glossa, Milton, B.S. in Ae.E.; Detail Draftsman, Boeing 
Aircraft Co. 
Harvey, Jacque W., B.S. in Ae.E.; Stress Analyst, Bell Air- 
craft Corp. 
Hickey, Whalen K., Draftsman, Airplane Div. (Buffalo Plant), 
Curtiss-Wright Corp. 
Hughes, George W., Jr., M.E.; Liaison Engineer, Consolidated 
Aircraft Corp. 
Jue, George, M.S. in Ae.E.; Stress Analyst, The Zap Corp. 
Kantrow, Harold, A.B.; Scientific Aide, Langley Memorial 
Aero. Lab., N.A.C.A. 
Lamar, William Edgar, B.S. in Ae.E.; Lt. U.S. Army Air 
Corps. 
Lewis, Richard B., B.S. in Ae.E.; Flight and Apprentice 
Engineer, Pan American-Grace Airways, Inc. 
Lorenz, Sigmund J., Drafting, The Glenn L. Martin Co. 
Mayer, Norman J., Aero. Engineer, Goodyear Aircraft Corp 
McMillan, William P., B.S. in Ae.E.; Student Instructor, 
1.8. Army Air Forces. 
Moore, Wesley Dean, Layout Draftsman, El Segundo Div., 
Douglas Aircraft Co. 
Mortlock, Walter E., B.S. in Ae.E.; Power Plant Analysis 
Engineer, Boeing Aircraft Co. 
Moses, Raymond E., Engineer, Airplane Div. (Columbus 
Plant), Curtiss-Wright Corp. 
Newmark, M. F., Jr., B.S. in Ae.E.; Inspector, Consolidated 
Aircraft Corp. 
Pellerin, Earl H., B.S. in-Ae.E.; Loftsman, Consolidated Air- 
craft Corp. 
Plommer, Wm. D., Design Engineer, Airport Div., Mackenzie 
Air Service Ltd. 
Podolsky, Arnold C., B.S. in Ae.E.; Stressman, Vought- 
Sikorsky Aircraft Div., United Aircraft Corp. 
Randall, Fred W., Jr., B.M.E.; Apprentice Engineer, Vought- 
Sikorsky Aircraft Div., United Aircraft Corp. 
Rechlin, Floyd, B.S. in Ae.E.; Aero. Engineer, Consolidated 
Aircraft Corp. 
Richards, Cyril E., B.S. in Ae.E.; Inspector, Consolidated Air- 
craft Corp. 
Soennichsen, Matthew E., B.S. in Ae.E.; Ensign, U.S.N.; 
Officer Asst. (Radio Elec. Lab.), U.S. Naval Aircraft Factory, 


— 


Philadelphia. 
Speer, Edmund Magnus, Jr., B.S. in Ae.E.; Test Inspector, 
Pratt & Whitney Aircraft Div., United Aircraft Corp. 
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Strickland, John T., M.S. in Ae.E.; Asst. Aero. Engineer, 
J.S. Army Air Corps. 
Swanson, Douglas. 
Tenedini, A. S., Student, Casey Jones School of Aeronautics. 
Vondenberg, Joseph L., B.S. in Ae.E.; Private, U.S. Army. 
Wall, Thomas F., Draftsman, Douglas Aircraft Co. 
Warshaw, Seymour, B.C.E.; Jr. Naval Architect, U.S. Navy 
Yard, Philadelphia. 
Winkler, Karl J., B.S. in Ae.E.; Jr. Engineer, Curtiss Pro- 
peller Div., Curtiss-Wright Corp. 
Wisely, Joseph A., B.S. in Ae.E.; Jr. Engineer, Curtiss Pro- 
peller Div., Curtiss-Wright Corp. 


NECROLOGY 
Hart O. BERG 


Hart O. Berg, Benefactor Member of the Institute, died in 
New York on December 9, 1941, after an illness of several weeks. 

Due to his early connection with the work of the Wright 
brothers and his continued interest in aviation, Mr. Berg had 
kept a large collection of documents on their work. In 1939 he 
presented this collection to the Aeronautical Archives of the 
Institute. They contain the most complete first hand record of 
the early flights of the Wright Brothers in Europe. 

Hart O. Berg was born in Philadelphia, on March 23, 1865. 
He received an engineering education in Liege, Belgium, and 
worked for several years in a Belgian gun plant. He accepted, 
in the early nineties, the position of Assistant Manager of the 
Colt Patent Fire Arms Company’s factory at Hartford, Conn. 
During his two years at the Colt plant the manufacture of the 
old Gatling gun was discontinued and the new Browning machine 
gun was developed under his direction. 

He then went back to Belgium where for six years he was suc- 
cessively Assistant Manager and General Manager of the Fab- 
rique Nationale d’Armes de Guerre at Herstal, Liege (Belgian 
National Gun Works), building Mauser rifles. He brought 
John Browning to Liege and manufactured the first Browning 
pistols at the Belgian plant before they were produced com- 
mercially in the United States. He also traveled extensively 
and obtained important gun and munition contracts in almost 
every country of the world for the Fabrique Nationale. 

In 1899 he became Chief Engineer for Clement-Bayard in 
Paris, and erected there the world’s first large automobile plant. 
He represented the United States as American Member of the 
Machinery Jury at the Paris Exposition in 1900. 

In 1902 he organized the Berg Automobile Company and 
built the ‘““Berg,’”’ one of the first large four cylinder cars con- 
structed in America. 

Shortly afterward he became associated with Captain Simon 
Lake, and obtained, early in 1904, an order from the Russian 
Government for six submarine boats. He spent three years in 
Russia building ten submarines for the Russian Navy, the 
last four being of a larger type than any previously built, each 
mounting for the first time a large calibre gun. He also built, 
with the late Lewis Nixon, ten torpedo Boats at Sebastopol. 
Afterward he built, with Simon Lake, two submarine boats in 
Pola for the Austrian Government. 

In 1908 Mr. Berg brought Wilbur Wright to Europe and hired 
a race course near Le Mans, France, where in September of that 


year Wright made his first public airplane flights and so started 
practical aviation. He then organized in France the world’s 
first aviation company. Shortly afterward Wilbur Wright won 
the Michelin Prize of 20,000 francs (then $4000) for directed 
flight over a line of toy balloons 100 meters high. He organized 
the first flying field at Pau and later a flying field outside of 
Rome, and the following year took Orville Wright to Germany 
where the first airplane flights in Germany took place at the 
Tempelhof, Berlin. During this period he also brought to 
France Lieutenant Reilly-Scott, a West Point graduate, who 
won the first Michelin Prize for dropping bombs at a target from 
an airplane. 

In 1910 he became interested in a new method of obtaining 
sulphur from waste gases in smelter operations, and shortly after- 
ward was instrumental in erecting the first oil cracking plant in 
Europe at the mouth of the River Thames to increase production 
of gasoline. 

Later during World War I he designed, equipped, and installed 
ten ice plants for the making and distributing of ice to Field 
Hospitals on the Eastern Front. 

Since the War he was employed in metallurgical operations 
principally concerning the production of metallic magnesium 
and its commercial introduction. 

He had made Paris his home for many years up to 1939 and 
had his main offices and staff for general engineering and con- 
sulting work there as well as branch offices in London and New 
York. He made more than 146 trips across the Atlantic, be- 
sides encircling the world several times. 

Mr. Berg was a member of the American Society of Mechani- 
cal Engineers, of which he was Honorary Vice-President and their 
representative at the Memorial Ceremonies in honor of the 
French scientist Le Chatelier, member of the Engineers Club, 
New York, one of the Founders in 1898 of the Aero Club de 
France, also of the Automobile Club de France, member of the 
American Club of Paris, one of the Founders of the St. Cloud 
Country Club, member of the Royal Automobile Club, London, 
and one of the first members of the American Chamber of Com- 
merce in Paris. 

He was decorated twice by the French Government; in 1901 
as Chevalier, and in 1909 as an Officer of the Legion of Honor. 


TuLuio Levi-CIvITa 


Dr. Tullio Levi-Civita, foreign Associate Fellow of the In- 
stitute, died in Rome, Italy, on January 1, 1942, at the age of 69. 

He was one of the small group of distinguished foreign scien- 
tists invited to be among the Founder Members of the Institute. 
Dr. Levi-Civita was regarded as Italy’s greatest mathematician 
and had done important work in the fields of hydrodynamics, 
theoretical dynamics and pure geometry. His contribution to 
the theory of absolute differential calculus was said to have helped 
make possible Einstein’s general theory of gravitational relativity. 
His concept of parallelism in higher spaces was regarded as show- 
ing the way to the development of the new differential geometry. 

Dr. Levi-Civita was educated at the University of Padua, 
from which he received a Doctorate in Mathematics in 1894. 
He had served as Professor of Rational Mechanics at the Uni- 
versity of Rome since 1919. He was a member of the Pontifical 
Academy of Sciences, the Associazione Aerotecnica Italiana and a 
foreign member of the British Royal Society. 











